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Abstract 

o" 

Sea quark distributions in the nucleon have naively been expected to be generated perturbatively by gluon 
splitting. In this case, there is no reason for the light quark and anti-quark sea distributions to be different. 
^ . No asymmetries in the strange or heavy quark sea distributions are predicted in the improved parton model. 
Q \ However, recent experiments have called these naive expectations into question. A violation of the Gottfried 
' sum rule has been measured in several experiments, suggesting that u < d m. the proton. Additionally, other 
^ , measurements, while not definitive, show that there may be an asymmetry in the strange and anti-strange 
\ quark sea distributions. These effects may require nonperturbative explanations. In this review we first discuss 
' the perturbative aspects of the sea quark distributions. We then describe the experiments that could point to 
: nonperturbative contributions to the nucleon sea. Current phenomenological models that could explain some 
00 , of these effects are reviewed. 
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1 Introduction 



o ■ 

O ' 

■ 

O ■ Elastic scattering of the proton with electrons revealed that, unlike the electron, the proton had a finite 
size Further experiments with large momentum transfers showed that, in addition to excitation of 
\ the proton into resonance states such as the A"*", the proton could break up into a large amount of debris. 
^ ■ Beyond the resonance region, it appeared as though the virtual photon probe was no longer interacting 
with the proton as a whole but with point-like constituents of the proton. Early measurements of the 
structure functions suggested that they were independent of the momentum transfer, Q [El- Out of these 
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data emerged the parton model of the nucleoli [Q. In this picture, the proton has three valence quarks, 
two up, or u, quarks and one down, or d quark. If these valence quarks, bound in the proton, were its 
only constituents, there would be nothing further to tell. However, there is much more to the story. 
When later experiments probed lower values of the fraction of the proton momentum carried by the 
partons, x, scaling was found to be violated [|]. Additionally, if the proton only contained valence 
quarks at all values of Q, each valence quark should carry a third of the total proton momentum. This 
also proved to be untrue. The average x of the partons was less than 1/3 and the parton density was 
correspondingly larger in this x < 1/3 region. The low charged x partons contributing to the proton 
structure functions are quark-anti-quark pairs, known as sea quarks. 

The sea quarks are perturbatively generated from gluon splitting, the gluons being the mediators of the 
strong force. Since the gluons are charge-neutral, they do not directly contribute to the proton structure 
functions. When no interaction occurs, the sea quarks and gluons serve to dress the bare valence quarks. 
However, an interaction with a high energy probe such as a virtual photon breaks the coherence of the 
proton wavefunction and hadronizes the sea quarks and gluons. The higher the energy of the probe, the 
smaller the value of x that can be studied. Currently the smallest values of x have been reached in ep 
collisions at HERA with a center of mass energy of 310 GeV. These data show a strong increase 
of the sea quark and gluon densities at low x. 

Inherently, the generation of sea quarks at low x and Q is a soft process and cannot be calculated in 
perturbative Quantum Chromodynamics, QCD. Typically a minimum perturbative momentum transfer 
is chosen at which the valence quarks, sea quarks, and gluons are assigned distributions for further 
evolution to higher Q^. Since these distributions are not calculable, certain assumptions are made 
regarding the sea quarks. Until recently, it was always assumed that the sea quarks were all symmetric- 
there was no difference between up, down and strange sea quarks and no difference between e.g. the up 
and anti-up quarks in the sea. In the last few years these assumptions have been seriously challenged. 
These challenges are the subject of this review. 



2 Generation of the perturbative sea 

2. 1 The improved parton model 

The parton model can be applied to the calculation of any cross section with large momentum transfer, 
referred to as a 'hard interaction'. At high energies, the hadron and parton masses are neglected com- 
pared to the hard scattering scale Q. Partons from the initial state hadrons participating in the hard 
scattering have momentum xp^ where x is the fraction of the parent momentum carried by the parton, 
< a; < 1, and p^ is the four- momentum of the initial hadron. In the hard scattering, = since the 
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hadron mass has been neglected. A hard cross section calculated in the parton model from collisions of 
hadrons A and B to make final states C and D is schematically 



dXidX2 -—aij{XiP, X2p')(t)i/A{Xl)(t)j/B{x2)Dk{zi)Df[Z2} 

Zl Z2 



Here aij{xip,X2p') is the lowest order partonic cross section for partons i and j to make the final state 
kl i.e. gg — *■ qq for quark jet production. The probability density for finding parton i in hadron h with 
momentum fraction x is (j)i/h{x). To produce a final state hadron, C, from parton k or /, the partonic 
cross section is convoluted with a fragmentation function, D^{z), that describes the hadronization of 
the parton. The functions (pi/hix) which underlie all perturbative QCD calculations are the primary 
focus of this review. 

The most transparent picture of the parton model is deep inelastic scattering of an electron with a 
hadron. In the hadron, partons exist in a virtual state with definite momenta. When the electron 
approaches the hadron, due to Lorentz contraction and time dilation, the virtual state is frozen during 
the time it takes for the electron to cross the hadron. When the lepton transfers the squared four- 
momentum to the hadron via the virtual photon, the distance of closest approach between the 
lepton and hadron is ~ l/Q- The ratio of the "area" of the momentum exchange to the geometric area 
of the hadron, {1/Q'^)/ttRI, measures the locality of the lepton probe. When is large, this ratio is 
small, implying that the parton density in the hadron is so dilute that only one parton interacts with 
the lepton. The hadronization time scale of the partons produced in the hard scattering is taken to be 
long compared to that of the hard scattering. Therefore, both the initial and final states are decoupled 
from the hard interaction 

Most of the existing information on the nucleon sea has been obtained from deep-inelastic scattering, 
DIS, and lepton pair production by the Drell-Yan process. The DIS and Drell-Yan kinematics are briefly 
introduced, following the review of Ref. ||^. While the discussion remains at leading order for clarity, 
higher order calculations are available, see e.g. A typical deep-inelastic scattering reaction is 

l{pi) + h{ph)^l'{pi')+X (2) 

where / is a lepton such as an electron, muon or neutrino, h is typically a proton or a nucleon, and X is the 
hadronic final state. A virtual vector boson such as a photon or a is exchanged in the intermediate 
state. DIS is an inclusive process so that the kinematics can be specified by the measurement of the 
final lepton four momentum, The momentum transfer is space-like, q^ = Pi — Py and — = . 
The momentum fraction x, the Bjorken scaling variable, is the ratio of the momentum transfer squared 
to the energy transferred from the lepton to the hadronic system in the scattering, 

X = = . (3) 

2ph ■ q 2m,, z/ 
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The second equality holds in the target rest frame where u = Ei — Eii is the energy transferred from 
the initial lepton. The ratio of the energy transfer to the total lepton energy in the target rest frame, 
y, is also useful, 

Ph-q V 

y = = • 4 

Ph ■ Pi El 

The mass of the hadronic final state, W, can be large when x is fixed and is large since 

= ml + ^{l-x) . (5) 

X 

Depending on the initial lepton and the squared momentum transfer, the exchanged boson can be a 
photon, 7, or one of the electroweak vector bosons, W~^, W~, and Z^. When a 7 or Z° is exchanged, 
there is no charge transferred by the virtual vector boson and the reaction is said to be a neutral current 
process. The exchange of a or a W~, on the other hand, is through a charged current. The following 
discussion applies to both neutral and charged current processes. 

The DIS cross section can be separated into leptonic and hadronic components so that 



dcT = , / ... L^^iPi, q)W,APh, q) (6) 



2s|p//| 47r2(g2 _ rriyY 

where = 7, W^^ or with mass my. The vector coupling constants are = e, Cw± = fi'/2-\/2, and 
Czo = g/2cos9w where g = e/ sinOw and 9w is the Weinberg angle. The boson propagator is assumed 
to be oc (g^ — my)~^. 

The lepton tensor is 

^'^'^(p/,g) = ^Tr[^,r^,(#,-^)r^j (7) 

where the couplings of the vector bosons to the leptons are r{^^± = 7^, r|^+^ = 7^(1 — 7^), r|^__ = 

7'^(1 + 7^), and F^o = 7'^(cy — c^aI^) where Cy and are the vector and axial vector couplings for 
leptons and quarks of flavor /. Since this review focuses on values of where the virtual photon 
contribution is much larger than that of the Z^, the Z^ will no longer be discussed. The lepton spin 
average factor of 1/2 in Eq. (|^) is absent in neutrino scattering. 

The hadron tensor is 

W,^ip,,q) = 27,'J2ih''{Ph)\jrmX){X\j^{0W{p,))5Hp, + q-px) (8) 

a,X 
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where the nucleoli spin, a, is averaged over and the final state, X, has been summed over. Symmetries 
allow the hadron tensor to be written in terms of structure functions describing lepton scattering with 
point particles so that 

W,.ip,, q) = - [g,. - Q') (9) 

mi 

where parity conservation in strong interactions implies = for virtual photons. In the Bjorken 
limit, u, Q"^ oo for fixed x, the structure functions Wi are functions of x alone to within logarithmic 
corrections in Q^, 

F^\x,Q') = W,{x,Q'), 

F^\x,Q') = —W,{x,Q'), (10) 

F^'^{x,Q^) = —W,{x,Q^) . 
rrih 

The DIS cross section can then be expressed in terms of x, y, and F^^ as 



dxdy 



N' 



IV 



2 / \ 



'-2xFrix, Q'')+(l-y- ^) F^ix, Q') 

fill 



+ Sv (y-^]xFP{x,Q') 



2 

where E is the lepton energy, ruh is the hadron mass, and 6v = for 7,-1 for W~, and +1 for W'^. 
The normalizations are 



where a is the electromagnetic coupling constant. The lepton-hadron cross section in Eq. ( [Tl| ) is equiv- 
alent to lepton scattering from point-like objects. These objects are the partons. 

In the parton model, only charged partonic constituents of the hadron, the quarks and anti-quarks, 
couple to the electroweak currents at leading order. The DIS cross section is then proportional to the 
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elastic quark- lepton scattering cross section multiplied by the probability of finding a quark or anti-quark 
of flavor / and momentum fraction C,, 



The structure functions can then be directly related to the probability densities by 

1 



^m(^) = E / Y^li (f ) (0/7/^(0 + 0J/.(e)) (15) 



J " 

F^^x) = E / diF^ ( I ) (0^/,(O + 0j/,(O) , (16) 



/ 

where F^^ are the parton-level structure functions calculated at leading order. 



Ih — > IX where / = or yU^. It can be shown that 2F2^ {x) = {x) = Q'j5{l - a;) and that, 
subsequently, the hadronic structure functions are 

2xF^\x) = Ftix) =j:Q>ihMx) + 07/. (x)) . (17) 



/ 

The simplified notation 4>u/p{x) = Up{x), (f)u/p{x) = Up{x), etc. is now introduced for the partons in the 
proton. In this case, Eq. (|T^ is typically written as 

2xFr{x) = Fr{x) = E Q'Afpix) + 7p{x)) (18) 

/ 

where f = u,d, s,c, ■ ■ ■. Note that f{x) is the parton probability density while xf{x) is typically referred 
to as the parton momentum distribution. The u and d quarks are already present in the proton since, 
in the constituent quark model, the proton is a uud state where the valence quarks carry the bulk of 
the quark momentum at large x. The sea quarks do not contribute to the baryon number of the hadron 
since they are always produced in a virtual quark- ant i- quark pair. Therefore, for every Up there is a Up 
(denoted Up to avoid confusion with the total Up quark density). The u and d quark densities include 
both valence and sea quark densities so that 

Up{x) = Mp(x) + Up{x) (19) 
dp{x) = dlix) + d'p{x) . (20) 

The number of valence quarks in the proton sums to three. 
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1 



/ 



dxUp{x) = 2 



(21) 







1 



/ 



dx dp{x) = 1 , 



(22) 







two up valence quarks and one down valence quark. 

In addition to the proton parton distributions, DIS on targets such as deuterium or other light nuclei 
measures the neutron distributions. Typically, one assumes that the per nucleon structure function is 
e.g. for a deuterium target 



where the neutron parton densities are related to those in the proton by charge symmetry so that 

Up{x)=dn{x) dp{x) = Unix) 



Charge symmetry is also used to establish the valence quark sum rules for the neutron, as in Eqs. (|2T|) and 
(|22|) . The structure functions of other baryons are unmeasured but could be inferred from relationships 
like those of Eq. (0). Parton distribution functions have also been determined with limited statistics 
for mesons. In general, the valence quark densities of the meson are assumed to be Jq dxq'lj{x) = 2 
where u and d are the valence quarks when M = tt^. In this case, m^- is a component of the pion sea. 

An exchanged photon treats all charged partons identically. Thus the electromagnetic structure function 
of the proton, Eq. ([TB|), is the sum over all parton flavors. However, in neutrino-induced interactions, to 
determine the momentum fraction of the proton carried by the interacting parton, the reactions must 
proceed by charge current scattering for the final-state electron or muon to be detected. Thus, the DIS 
reactions are e.g. z/g/i e~X with an exchanged and Vgh — > e~^X with an exchanged W~. Now 
not all quarks and anti-quarks couple to the and W~. In the reaction Ueh — > e~X, the incoming 
Ue "decays" to W+ by i/^ ^ W+e'. Then ly+rf u, W+s c, W+u d, and W+c ^ s so that the 
couples to the d, s, u, and c quarks in the proton. Similarly, the W~ emitted in the —>■ W~e^ 
conversion can couple to the d, s, m, and c quarks in the proton. Therefore, for a neutrino beam with 
exchanged , the structure functions are 



Ff{x) = \{Fr{x)+Fr{x)) 



(23) 



Up{x) = dn{x) dp{x) = Un{x) 



(24) 



F2 ^(x) = 2x [dp{x) + Sp{x) + bp{x) + Up{x) + Cp{x)] 
FY'^^ix) = 2 [dp{x) + Sp{x) + bp{x) - Up{x) - Cp(x)] 



(25) 
(26) 
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Note that the b quark has been included in the proton sea while the still heavier top quark has not. 
Likewise, the structure functions determined by an anti-neutrino beam with exchanged W~ are 



-,W-p/ 



X] 



2x 



dp{x) + Sp{x) + bp{x) + Up{x) + Cp{x) 
dp{x) + Sp{x) + bp{x) - Up{x) - Cp{x) 



(27) 
(28) 



Some DIS data, particularly with neutrino beams, have been taken on nuclear targets heavier than 
deuterium. Therefore it is useful to write the electromagnetic, neutrino and anti-neutrino structure 
functions in terms of an isoscalar nuclear target. (In an isoscalar target, Z = N = A/2, where A is 
the nuclear mass number, Z is the number of protons and N is the number of neutrons.) Assuming no 
modifications of the parton densities in the nucleus, the resulting structure functions per nucleon are: 



x) = —x[up{x) + Up{x) + dp{x) + dp{x) + -(c(x) + c{x)) 

-Lo 
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(x) + Sp 


{x) + bp 


[x) + bp 


»)] , 












pW+N, 


[x) 


= x[up 


[x) 


+ Up{x) 


+ dp{x) 


+ dp{x) 


+ 2sp 


[x) 


+ 2bp 
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+ 2cp 


»] , 






= x[up 
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+ dp{x) 


- Up{x) 


- dp{x) 


+ 2sp 


[x) 


+ 2bp 




— '2cp 


W], 


pW-No 




= x[up 
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+ dp{x) 


+ dp{x) 


+ 2sp 


[x) 
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+ 2cp 


»] , 






= x[up 
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+ dp{x) 


- Up{x) 


dip (^^^ 


— 2sp 


[x) 


- 25p 


[x] 


+ 2cp 


[x)]. 



(29) 

(30) 
(31) 
(32) 
(33) 



Since measurements of the lepton pairs produced in the Drell-Yan process play an important role in 
the determination of the nucleon sea, the Drell-Yan process 0] is now introduced. At leading order a 
quark from the projectile annihilates with its corresponding anti-quark in the target, producing a virtual 
photon which then decays to a lepton pair with invariant mass M, gg — > 7* — > The Drell-Yan 

cross section involves a convolution of the parton distributions of both the projectile and the target 
hadrons. The partonic cross section for Drell-Yan production is 

-Q)5{s-M^) (34) 



dM 9M 

where s = X1X2S, Xi is the fractional momentum carried by the projectile parton, X2 is the fractional mo- 
mentum of the target parton, and s is the square of the hadron-hadron center of mass four-momentum. 
To obtain the hadroproduction cross section as a function of pair mass the partonic cross section is 
convoluted with the quark and anti-quark densities evaluated at scale M, 



1 

vra' 



dM 9M 



I dx,dx2 6{s - M2) ^ Q}[fp{xu M')f^{x2, 
/ 
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+ /p(xi,M2)/p(x2,M2)] 



(35) 



Besides deep-inelastic scattering and Drell-Yan pair production, prompt photon, heavy quark, quarko- 
nium, and jet production are among the processes calculable in perturbative QCD. To obtain the 
hadronic cross sections, the partonic cross section must convoluted with the projectile and target par- 
ton densities. For such a convolution to be reliable, the hadronic cross sections must be factorizable into 
the process specific short-distance cross section and the parton distributions which cannot be obtained 
from a first principles calculation and must be assumed to be universal |ll|], as in Eq. (|l|). Since the 
parton distributions are properties of the hadron wavefunction and thus nonperturbative in nature, a 
universal set of parton densities for a given hadron is typically obtained from fits to DIS and other data. 
How this is generally done is described in the next sections. 



2.2 Evolution 



Deep inelastic scattering data are available over a wide range of x and In the naive parton 

model, scaling was assumed, i. e. the parton distributions were thought to be independent of the of 
the virtual vector boson p[. However, further measurements 0] soon showed that this was not the case. 
To obtain a reliable set of parton densities at all x and Q^, a description of their evolution in x and 
must be available. While this is relatively simple at fixed x when is large, the small x behavior of 
the parton densities is more challenging, as will be discussed later. 

Only tree-level (leading order) calculations have been described so far. However, many processes have 
now been calculated at least to next-to-leading order in QCD. Next-to-leading order calculations involve 
both real and virtual corrections to the leading order diagrams. At any order in perturbation theory, 
the hadronic cross section factorizes into a perturbatively calculable partonic cross section and the 
parton distributions in the hadron, as in Eq. (|I]). Factorization provides a scheme for separating the 
short distance partonic cross section from the nonperturbative parton distribution functions, generally 
referred to as a factorization scheme. The momentum scale at which the separation occurs is the 
factorization scale. 

There are two commonly used factorization schemes, the DIS and the MS schemes. In the DIS scheme, 
FJ^'^ is defined so that all corrections are absorbed into the quark and anti-quark distributions order- 
by-order in perturbation theory. Thus F^^ is trivial. However, only F^^ enjoys this simplicity while all 
other structure functions must be corrected order- by-order in a^. Additionally, the gluon distribution is 
not fixed in this scheme. The more commonly used MS or minimal subtraction scheme defines the parton 
distributions in terms of parton creation and annihilation operators acting on the hadron wavefunction. 
The next-to-leading order corrections to the hard scattering are somewhat more complex but the MS 
scheme is more physical and includes the gluon through evolution. 
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In addition to the factorization scheme and scale dependence, another scale, the renormalization scale, is 
required to remove the ultraviolet divergences that occur in the calculation of the partonic cross section. 
In the MS scheme, a uniform mass or energy scale is assigned to all diagrams required to calculate the 
cross section, helping to regulate the divergences. Such a scale assignment automatically introduces a 
running coupling Q!s(//|j) which allows the cross section to remain independent of the renormalization 
scale fiji- 

There is no reason why the factorization and renormalization scales should be different and they are 
indeed usually chosen to be the same. Choosing the two scales to be equal in e.g. the MS scheme ensures 
that the parton distributions can be calculated with the local operators and divergences in the cross 
sections can be handled simultaneously. The cross section, a physical quantity, should not depend on 
the choice of scheme or scale and would not if it could be calculated to all orders in perturbation theory. 
However, since the expansion of the cross section in powers of the strong coupling constant is truncated 
at finite order, a dependence on the scale and scheme remains. The scale at which the hard scattering 
occurs is then generally the one at which the parton distributions are evaluated. 

An important byproduct of factorization is that masuring an observable such as at one value of 
such as Qq enables one to predict the behavior of at another value of via evolution. In the 
evolution from Ql to Q^, a parton of momentum fraction x, probed by a virtual photon with momentum 
transfer Ql, can split into two partons with lower momentum fractions, y and z, with y + z = x. For 
a virtual photon to resolve a charged parton of momentum y, y < x, it must have a larger momentum 
transfer than Qq, > Q^. When x and are not too small, the Dokshitzer-Gribov-Lipatov-Altarelli- 



Parisi (DGLAP) equations [12| describe how the parton distributions are determined at > Qq 



starting from an initial distribution at a low scale Q^. Evolution in links the charged quark and 
anti-quark distributions with the neutral partons, or gluons, since quarks emit gluons while gluons can 
split into qq pairs or into two gluons. The gluon is thus defined through its evolution in the MS scheme. 
At leading order, the DGLAP equations are 



cilogQ2 



dqi{x,Q'^) _ as{Q^) f dy 
dlogQ^ 

dg{x, Q2) 



2vr { y 
as{Q^) } dy 



2v: J y 



(36) 
(37) 



where the sum i = 1 . . .2nf runs over quarks and anti-quarks of all flavors. The probability that parton 
i has emitted parton j is proportional to the splitting functions Pij. The strong coupling constant, a^, 
appears because of the quark-gluon and gluon-gluon vertices. The parton distributions in Eqs. ( [I^ ) - 
(|28|) all have had the Q"^ dependence suppressed. The suppression of the argument of the parton 
distribution functions is common throughout this review. 



While the evolution of the parton distributions can be described by perturbative QCD, some initial set 
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of parton distributions must be defined at scale Q^. This initial scale is somewhat arbitrary but can 
be thought of as the boundary between perturbative and nonperturbative physics. The initial scale is 
usually in the range 0.5 < Qq < 5 GeV^. 



2.3 Experimental studies 



Although a great deal of precise information is generally obtainable from DIS, particularly since the 
parton densities of only one hadron are involved, DIS does not distinguish between quark flavors. Other 
processes are also needed to identify the effects of individual parton flavors. A large collection of data 
are used to fit the parton distributions at the initial scale Qo- The DGLAP equations then evolve the 
parton densities to higher scales. A list of some relevant measurements and which parton distributions 
they illuminate is given below: 

- F2^{x,Q'^) from e^p interactions at the HERA collider Pfl, charged partons at small x; 

- F^^x.Q'^) and F:f^{x,Q'^) from ^l beams by NMC [H and BCDMS H, charged partons; 



F^^^ {x^Q"^) and xFY^^^ {x,Q'^) from v beams on nuclear targets pr5|JT6|JT7[] , charged partons; 



2 y-^i^ ) "■^'■^ -^^ 3 

- prompt photons, qg — > 7g, gluon distribution; 

- Drell-Yan production, qq ^ ^ , sea quark distributions; 

- and Z° production in pp colliders, u and d distributions — e.g. awi/o'zo oc d/u; 

- Z° pole measurements in e~^e~ collisions |TB[, intrinsic scale, Aqcd, of 



Note that combinations of F^^ and F^^ can pick out certain parton densities. For example, using 
Eqs. (HD-dlD, xF^^'' + xFY~^ + F^*"^ - F^'" = xd;{x,Q^). In fact, if both p and n targets are 



used with the assumption of charge symmetry and an SU(3) flavor symmetric sea, Up = dp = Sp with 
Cp = bp = 0, the parton distributions are overdetermined and consistency checks can be made. It turns 
out, as discussed in this review, that there is no sea quark SU(3) flavor symmetry. 

Some characteristic data on F2^ from muon beams at 120, 200, 240, and 280 GeV on hydrogen targets 
at CERN and the neutrino structure functions F^ ^ and xF^ ^ from u and 17 beams on an iron 
target at Fermilab are shown in Figs. |l| and ^ respectively. Several values of x are selected and the 
dependence of each x bin is shown. Note that in the intermediate x range, 0.125 < a; < 0.35, little 
dependence is observed, prompting the scaling assumption from the early data taken in this x region 
0. More recently, the HERA ep collider at y/s = 310 GeV has provided data over a wider range of x 
and than previously available. A sample of the HERA data is shown in Fig. ^. The scaling violations 
only somewhat apparent in the fixed target experiments are more obvious at HERA. 

Since pure neutrino beams are unavailable, neutrino experiments with secondary beams collect data 
from both neutrino and anti-neutrino-induced events from all possible neutrino energies. The resulting 
measured structure functions are a weighted average of the u and u induced events, e.g. F^'^^° and 
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Fig. 1. The proton structure function as a function of for muon beams with energies 120, 200, 240, and 
280 GeV measured by the European Muon Collaboration ||T^ . From top to bottom the values of the x bins are: 



0.0175, 0.03, 0.05, 0.08, 0.125, 0.175, 0.25, 0.35, 0.45, 0.55, 0.65, and 0.75. To separate the data, F:]^ in each x 
bin is scaled by a factor of 1.5 from the next higher x bin. Therefore, only the highest x data (lowest points) 
have the correct scale. Only the statistical uncertainty is shown. 

pW No j^p^ Ny/[Ny + Ny) Is thc fraction of DIS events induced by neutrinos, the experimentally 
determined structure function is 



:/3Ff^^»(x) + (l 



X) 



[Ff''^«(x) + F2"(x)] 



W+No 



XI 



X)] 



(3J 



If /3 = 1/2 or F^ 



W+No 



^2 



implying Sp{ ), the second term in Eq. ( PBD vanishes. From the 



CCFR data ||20|, j3 = 0.83 |^ so that F2 ° is dominated by F2 ° or neutrino-induced production. In 
an analysis of opposite sign dilepton production, they found that to next-to-leading order, Sp{x) = Sp{x) 
within their experimental errors (These results will be discussed in more detail in Sec. 3.2.1.) 
However, since (3 > 1/2, the associated uncertainties are significant. By comparing the neutrino and 
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Fig. 2. The nucleon structure functions -Pj^ ^ (a) and xF^ ^ (b) as a function of from neutrino-iron 

deep-inelastic scattering. The x bins from top to bottom are: 0.015, 0.045, 0.08, 0.125, 0.175, 0.225, 0.275, 0.35, 
0.45, 0.55, and 0.65. To separate the data, the structure functions in each x bin are scaled by a factor of 2 
from the next higher x bin. Therefore, only the highest x data (lowest points) have the correct scale. Only the 
statistical uncertainty is shown. 

electromagnetic structure functions, one can form the relation. 



X] 



x] 



-X[Sp[X) + Sp[X)\ 

+ ^(2/3-l)Mx)-Sp(x)] 



(39) 



If Sp{x 



SpyX J 



then Eq. (p9D could be used to determine the strange quark distribution. However, 
the strange quark distribution extracted in this fashion does not agree with independent experimental 
results 1^ . 



To obtain the parton densities, a global analysis of all available data is performed. This analysis involves 
taking data from many different processes measured over a range of x and and making a multi- 
parameter fit to all the data, as described below. First, it is necessary to solve the coupled integro- 
differential evolution equations, Eqs. (BBf) and (p^) and at leading order or their next-to-leading order. 
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Fig. 3. The proton structure function FiJ^ as a function of at the HERA colhder (26.7 GeV electrons on 
820 GeV protons) §. From top to bottom the values of the x bins are: 0.000178, 0.000261, 0.000383, 0.000562, 
0.00075, 0.000825, 0.00133, 0.00237, 0.00421, 0.0075, 0.0133, 0.0237, 0.0421, 0.075, and 0.133. To separate the 
data, the structure function in each x bin is scaled by a factor of 2 from the next higher x bin. Therefore, only 
the highest x data (lowest points) have the correct scale. Only the statistical uncertainty is shown. 

NLO, counterparts, numerically to be able to calculate the parton densities at the appropriate to the 
measurement. The data to be used in the fit are then chosen. Typically only data sets that can provide 
the best constraints on the parton densities are included except when the data indicate significant new 
physics at work. Since the evaluation is done at fixed order in perturbation theory, a factorization scheme 
is chosen and all the evaluations of the parton densities are done for the same order of perturbation 
theory. It is very important that this be done consistently when calculating all processes considered. A 
starting distribution at scale Ql, typically of the form 



xf{x,Ql) = Aox^Hl-x)^2P(x) 



(40) 



is assumed where P{x) is a smooth polynomial function, the exponent Ai determines the small x 
behavior, and A2 governs the large x behavior. These initial distributions are then evolved to the 
appropriate scale of the data and used to calculate the cross section or structure function at the same x 
and as the data. The goodness of the fit, its x^? is calculated from a comparison to all the data in the 
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fit and the parameters in e.g. Eq. (^) adjusted in each iteration to minimize the x^- The final fits are 
generally made available either in tables with some interpolation formula used to obtain the densities 
at any x and or as parameterized functions of x and Q^. Parton distribution functions determined 
in this manner are only as good as the data they fit and may fail to match new data outside the x and 
range of the old measurements. As new data are taken, the parton distribution functions are thus 
continually updated. 

An example of updating the parton distribution functions as new data appear is the changing starting 
assumption regarding the degree of fiavor symmetry in the sea. In early global analyses, an SU(3) fiavor 
symmetric sea was assumed, see e.g. Ref. [^. However, it was shown that the strange sea is smaller 
than the light quark sea so that 2sp/ {up + dp) ~ 0.5 |1T7|J2^ . Thus, the sea could at most be SU(2) fiavor 
symmetric. It has been further shown that the SU(2) fiavor symmetry may not hold, as discussed in 



the next section. Newer fits of the parton distribution functions such as those in Refs. ||25|j26|p7| include 
these results in the sea quark analysis. 

Several different groups have been involved in the global analysis of the parton distributions. The CTEQ 
(Coordinated Theoretical-Experimental Project on QCD) collaboration |^5|] and Martin, Roberts, and 
Sterling (MRS) and their collaborators have produced a variety of parton distribution functions 



which improve with each round of data. Gliick, Reya, and Vogt (GRV) |]2^ have taken a rather different 
route, beginning at a low initial scale with valence-like sea quark and gluon distributions without trying 
to make a comprehensive fit. A library of all available parton distribution functions can be found in 
the program package PDFLIB available from the CERN program libraries Care must be taken 



however to use only the most recent sets, matching the appropriate order of the parton densities and 
hard scattering cross section in scale and scheme. 



In Figs. |, ^ and|^, the parton distributions obtained in recent global fits to NLO by CTEQ ||2^, MRST 
and GRV are shown as a function of x for two different values of Q^. In each case, the sea 



quark and gluon evolution is stronger than the valence evolution. The evolution is faster at low x than 
at intermediate and large x. While the Up and dp distributions are different for x > 0.01, as discussed 
shortly, they are nearly equal at low x. The strange quark distributions are generally smaller than the 
light quark distributions but by x ~ 10""^, Sp ~ Up. The gluon distribution, least constrained by data, 
evolves more rapidly in than do the charged parton distributions. 



2.4 Aspects of the small x sea 



As already discussed, the parton distributions are determined in practice from fits to DIS and related 
data by parameterizing the distributions at a Qo large enough for perturbative QCD to be applicable. 
Then the distributions are evolved to higher using the DGLAP equations, Eqs. (P^)-(P^). Evolution 
to higher when x is not too small should give reliable results. However, extrapolation to small x 
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Fig. 4. The CTEQ5 NLO MS scheme proton parton distributions are given at = 4 and 100 GeV^. The up 
and down valence distributions are shown in the solid and dashed lines respectively in (a) and (b). The up, 
down and strange sea distributions and the gluon distributions are given in the solid, dashed, dot-dashed and 
dotted curves respectively in (c) and (d). 



is more difficult. Early expectations based on Regge theory suggested that the sea quark and gluon 
distributions should become constant as x — 0. The HERA data have clearly shown that this is not 
true, as seen in Fig. ^. In this section, the way this behavior can arise is explained. 



Gluon ladder diagrams such as the one in Fig. taken from Ref. |2^, represent the behavior of 
DIS scattering. In the regime of finite x and large Q^, the DGLAP equations hold and the dominant 
contribution to the evolution of the sea and gluon distributions come from the region of strongly ordered 



transverse momentum, Q" 



When X is small and remains large, strong ordering in the longitudinal momentum, a; <^ ^„ <^ ■ ■ ■ <^ 
.^1, is as important as the transverse momentum ordering. The probability of emitting the i^^ gluon on 
the ladder is proportional to as[dxi/xi][dkfrp/kiT]'^. The nested integrations over x and kr generate the 
leading order behavior of the sea quark and gluon distribution, assuming is constant. 
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Fig. 5. The MRST NLO MS scheme proton parton distributions are given at = 4 and 100 GeV^. The up 
and down valence distributions are shown in the soUd and dashed Unes respectively in (a) and (b). The up, 
down and strange sea distributions and the gluon distributions are given in the solid, dashed, dot-dashed and 
dotted curves respectively in (c) and (d). 



xf{x, Q"^) ~ exp 



L TT \X, 



1/2 



(41) 



known as the double leading logarithm approximation, DLLA. The name DLLA implies that only the 
terms proportional to ln(l/x) InQ^ multiplying each power of as are kept in the expansion of Eq. (pj). 
In this case, xf{x, Q"^) increases as x — 0. The increase is faster than any power of ln(l/x) but slower 
than any power of x. 

In the case of small x and finite Q^, only the longitudinal momentum is now strongly ordered and the 
full kx phase space must be integrated over rather than retaining only the leading terms in InQ^. A 



counterpoint to the DGLAP equations at small x is the Lipatov or BFKL |^ equation that evolves a 
starting distribution /(xq, k"^) at Xq downwards in x by integrating over the kx phase space. 



X 



dfjx, kl) 
dx 



dk'SK{kT,k'T)f{x,k'S) . 



(42) 
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Fig. 6. The GRV 94 HO NLO MS scheme proton parton distributions are given at = 4 and 100 GeV^. The 
up and down valence distributions are shown in the sohd and dashed hnes respectively in (a) and (b). The up, 
down and strange sea distributions and the gluon distributions are given in the solid, dashed, dot-dashed and 
dotted curves respectively in (c) and (d). 

The kernel K includes real and virtual gluon emission. The Lipatov equation is related to the sea 
quark or gluon distribution before the integration over transverse momentum has been performed, 
f{x, k^) = d{xg{x, Q^))/^ In 1^2=^2 . The hadronic cross section is typically the partonic cross section 
convoluted with an x-dependent parton distribution, as in Eq. (|l|). However, in this case, the partonic 
cross section is convoluted with a fcr-dependent parton density so that cr ~ / cix [dk'^/k^]f{x, k^)a. If 
as is fixed, an approximate analytical solution may be found for x —>■ 0, 



xf{x,Q^) h{Q^)x-^ ; A 



12a., 



In2 . 



vr 



(43) 



The predicted value of A is ~ 0.5 when is fixed. Allowing to run with introduces a dependence 
on the cutoff of the k^ integration. 



The increase in the sea quark and gluon distributions at small x in Eq. (^3]) cannot continue indefinitely 
because the hadron will become so dense that the partons can no longer be considered to be free. The 



growth must eventually be suppressed, most likely by gluon recombination UST]]. It is unknown at what 
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Fig. 7. Ladder diagram corresponding to evolution of parton distribution functions. At large Q^, DGLAP 
evolution sums diagrams with strong ordering in kx, ^ kj, ^ k'^rp ^ • • • 3> kfj, ^> Qo- small x, 
the leading contribution comes from ordering in ^, x ^ ^ • • • ^ but with the gluon kx unordered. 
Reproduced from Ref. [29| with permission of the publisher. 



value of X this "shadowing" becomes important. For example, an anti-quark with distribution q{x, Q^) 
has an apparent transverse size of 1/Q but a smaller longitudinal size 1/xp with proton momentum p 
in a frame where xp 3> Q. The scattering cross section of the sea quarks and gluons is proportional to 
the square of the transverse size, a ~ as{Q'^) / Q'^ ■ When the area covered by the total number of sea 
quarks and gluons, n, scattering with the cross section a, approaches the transverse area of the proton, 
vrrp, or na 
equations p^ , |33| . 



~ vrr^, shadowing is no longer negligible and introduces an additional term in the evolution 



In practice, the parton distributions are only fit with DGLAP evolution with Ai chosen in accord with 
the BFKL behavior. The starting sea quark and gluon distributions in Eq. (|iO| ) have been found to have 
Ai ~ —0.3. While this is a faster growth than Regge theory suggested, it is lower than the idealized 
BFKL result of Ai ~ —0.5. It is difficult to see from these results whether or not shadowing is important 
because a full next-to-leading order evolution scheme incorporating DGLAP, DLLA, and BFKL is not 
available over the entire x and range. 
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2.5 Perturbative sea summary 

The naive parton model and its extensions have been extremely successful in describing a large number 
of processes. Both the data and the NLO perturbative calculations used to model the data are now 
extremely precise over a wide range of x, from 2 x 10"^ to 0.75. The parton distribution functions 
extracted from global analyses by a number of groups agree well with the data and with each other over 
their range of validity and will continue to be refined as new data are taken. 



3 Possible nonperturbative contributions to the sea 

So far, only the generation of the perturbative sea has been discussed. However, since the parton 
distributions are fit rather than derived from first principles and the perturbative evolution cannot go 
backwards and remain stable, the true nature of the nucleon sea is, at its heart, unknown. While the 
valence and sea quark distributions have been measured as a whole in DIS, and the valence and sea 
quark distributions can be separated from one another, the individual sea quark distributions are more 
difficult to tease out from bulk measurements. Therefore, simplifying assumptions have been made 
about the proton sea. Recently, experiments have been challenging the standard assumptions of the 
parton model such as SU(2) flavor and charge symmetry. The rest of this review is devoted to some of 
these results and their interpretations. 



3.1 The Gottfried sum rule^ vs. d 

Assuming flavor symmetry, proton and neutron parton distributions can be related. Valence quark 
counting suggests that the neutron d distribution and the proton u distribution should be equal, dn = Up 
at all X and Q^. Similar assumptions can be made for the other parton distributions, see Eq. (0). 
Therefore, unless flavor symmetry breaking is discussed, u, d, ■ ■ ■ can be assumed to refer to the proton 
parton distributions. The difference between the proton and neutron structure functions can then be 
written as 

i (Frix) - F^ix)) = 1 K(x) - d^ix)] + ^ [u{x) - d{x)] . (44) 
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Using Eqs. (^) and ( ^21) for the proton valence quark distributions, after integrating both sides of 
Eq. (HI) over x, one obtains 



j'^[Fr{x)-Fr{x)] = - + -jdx[u{x)-d{x)] . (45) 



If the sea is flavor symmetric, u = d and the second term vanishes. This is the Gottfried sum rule |p4|| : 
l^[Fn^)-Frix)] = l . (46) 



The assumption of light anti-quark flavor symmetry is a serious one and if the sum rule is found to be 
violated, the nucleon sea could be SU(2) flavor asymmetric. 



3.1.1 Current experimental studies 

Because the integral of Eq. (^) is over all x, the small x region could have a signiflcant contribution 
to the sum rule. Therefore the determination of any possible violation is only partial and depends 
strongly on the extrapolation of the parton distributions to low x. This was especially true when the 
proton parton distribution functions were not well constrained at low x by other measurements. The 
minimum x is restricted by the lepton-beam energy, E, to be Xmin = Q'^/2ME, where should not 
be smaller than a few GeV^ in order for the parton model to be applicable. To quantify the discussion, 
the Gottfried integral can be deflned as 

/G(x„.in,w)= / —[Fr{x)-Fr{x)] . (47) 

J X 

Generally two values of Iq are given for each measurement, one with Xmin and Xmax limited by the 
experiment, the other extrapolated over the entire x range. The integral over all x depends on how the 
missing small x region is accounted for. The results are summarized in Table |l]. When the Gottfried 
integral is given more than once for the same collaboration, the second result is obtained from a 
reanalysis of the earlier data. 

Typical DIS targets used to measure the structure functions are hydrogen and deuterium. The neutron 
distribution is then generally extracted taking into account the momentum smearing effects due to 
Fermi motion in the deuteron. However, the NMC analysis obtained the difference from the relation 

^2 -^2 -^^2 YTWTW ' 
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Measurement 



mm: •^max 



/g(0,1) 



sLAc m 



EMC |3§ 
EMC 1^ 
BCDMS [| 
NMC H] 
NMC M] 



0.02 0.82 0.20 ±0.04 

0.03 0.65 0.18 ±0.01 ±0.07 

0.02 0.80 0.197 ±0.011 ±0.083 

0.06 0.80 0.197 ±0.006 ±0.036 

0.004 0.80 0.227 ± 0.007 ± 0.014 

0.004 0.80 0.221 ±0.008 ±0.019 



0.24 ±0.02 ±0.13 



0.235 



-0.110 
-0.099 



0.240 ± 0.016 
0.235 ± 0.026 



Table 1 

Deep inelastic scattering measurements of the Gottfried sum rule. The Gottfried integral, Eq. (^), is given for 
each measurement and for the full x range when available. 



where the ratio Fa'^'/F^^ = 2F2^ / ~ 1 is determined from the NMC data [|38| , |39| . The absolute value 
of is fit from data. Nuclear corrections such as Fermi motion were not taken into account. The 
difference measured by NMC is shown in Fig. 

The data in Table |I] indicate /g(0, 1) < 1/3, a clear violation of the Gottfried sum rule. When the first 
data appeared, the low x parton distribution functions were not known from other sources. Therefore 
attempts were made to satisfy SU(2) flavor symmetry by changing the shape of the valence u and d 
distributions in the unmeasured region, see e.g. Ref. ESl. However, the lower x NMC data showed 



that this idea was untenable. Other attempts to understand the sum rule violation centered around the 
extraction of the neutron distribution from deuterium or other light nuclear data. The nuclear parton 



in nuclear targets, known as shadowing, which was typically neglected in deuterium. The effect has 
been studied in a variety of models including vector meson dominance, Pomeron and meson exchange 



distribution functions are different from free proton distributions [44|. At small x, there is a depletion of 



models and has been found to contribute only 5/g(0, 1) = —0.02 to the Gottfried sum rule. Note 

that since the shadowing contribution, Slc{0, 1) above, is negative, correcting for the effect in deuterium 
increases the Gottfried sum rule violation. Therefore, an SU(2) flavor asymmetry in the nucleon sea is 
the most likely explanation of the observed effect with d > u. 

Unfortunately DIS data can only address the difference between the light anti-quark distributions, 
usually written a.s d — u for a positive quantity. The individual anti-quark distributions cannot be 
studied separately. An independent measurement of d/u is needed. The Drell-Yan process provides an 
elegant way to obtain the ratio d/u by comparing Drell-Yan production on proton and deuterium targets 
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Fig. 8. The difference F2 — F2 measured by NMC |39[] at Q = 4 GeV compared witli calculations using 



the MRS (R2) g (dashed) and the MRST [||] (sohd) parton distributions. The E866 Q Drell-Yan results 
on the sea quark contribution to the difference are shown at the same Q^. For each set of parton distributions, 
the upper curve is the valence contribution to the difference, the lower curve is the sea contribution to the 
difference and the center curve is the sum of the two. Reproduced from Ref. [^] with permission from the 
Annual Review of Nuclear and Particle Science, Volume 49 @ 1999 by Annual Reviews. 



+ (4m(xi) + d{xij){u{x2) + d{x2)) 

assuming charge symmetry. When Xi is large, Il(a;i), d{xi) ~ so that 
1 



(50) 



a D ^ -(4m(xi) + d{xi)){u{x2) + d{x2)) 
y 

cTpp ^ ^(4'u(xi)m(x2) + d{xi)d{x2)) . 
The ratio of the pD and pp Drell-Yan cross sections when xi ^ X2 is 



(51) 
(52) 



DY 
'pD 



2aDY 

^ pp 



l + l[d{x,)/uix,)] / ^ dix2) \ 



2 1 + i [(i(xi)/u(a;i)] d{x2)/u{x2) 



u{x2) J 



(53) 
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Table 2 

Ratio of Drell-Yan cross sections in deuterium to hydrogen, (T^/2ct^ and the ratio d/u as a function of X2. 
The average of the kinematic variables for each bin is also tabulated. The units of {pt) and {M^+^-) are GeV. 
Note that {xi) = {xp) ± {x2)- Only statistical errors are shown. There is a 1% systematic uncertainty in all 
points. Modified from Ref. [^]. Copyright 1998 by the American Physical Society. 

When d = u, the ratio is 1. 

The first test of SU(2) flavor symmetry with the Drell-Yan process was made by the Fermilab E772 
collaboration. They compared the production of Drell-Yan muon pairs from isoscalar targets to that 
from a neutron-rich target and set constraints on the difference between u and d in the range 0.04 < 
X < 0.27 [0. Later, the CERN NA51 collaboration carried out a comparison of the Drell-Yan yield 



from hydrogen and deuterium at a single value of x with a 450 GeV proton beam and found 

= 0.51 ±0.04 ±0.05. (54) 



Up 

dp 



(a;>=0.18 



The most recent test has been made by the Fermilab E866 collaboration. They measured the Drell-Yan 
yield from an 800 GeV proton beam on liquid deuterium and hydrogen targets and extracted d/u and 
d — It in the proton for 0.020 < x < 0.345. The resulting ratio of the Drell-Yan cross section per nucleon, 
Eq. (0), is given in Table ^ as a function of Xi and X2 [ pn| , |5T| . The J/ip and T resonance regions, 
< 4.5 GeV and 9.0 < M^+^~ < 10.7 GeV respectively were excluded from the analysis. The 
data show that the Drell-Yan cross section per nucleon in pD interactions exceeds the Drell-Yan pp 
cross section over a range of X2- 
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Since all the data do not satisfy xi ^ X2, as necessary for Eq. (|53D to hold, the ratio d{x2)/u{x2) 
was extracted iteratively by calculating the leading order Drell-Yan cross sections using a set of parton 
distribution functions as input and adjusting d/u until the calculated ratio agreed with the measured 
one. In this analysis the light anti-quark sum, d + u, the valence, and heavy quark distributions obtained 
in global analyses of CTEQ 4M |^ and MRS (R2) |^ were assumed to be valid. The extracted d/u 
ratio is shown in Fig. along with the corresponding CTEQ 4M and MRS (R2) ratios. For consistency. 
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Fig. 9. The ratio of d/u in the proton as a function of x extracted from the Fermilab E866 cross section ratio 



pi] ] along with the NA51 result |49]. The solid curve shows the CTEQ 4M ratio and the dashed curve is the 
MRS (R2) ratio. The MRST ratio is shown in the dotted curve. The E866 error bars are statistical only with 



an additional systematic uncertainty of ±0.032 not indicated. Reproduced from Ref. |42| with permission from 
the Annual Review of Nuclear and Particle Science, Volume 49 (c) 1999 by Annual Reviews. 

when xi < 0.345, the projectile ratio d{xi) /u{xi) was assumed to be equal to that in the target proton 
while for xi > 0.345, d{xi)/u{xi) = 1 was assumed. When this high xi assumption was relaxed, there 
was a negligible difference in the low X2 results and only a 3% change in the largest X2 bin. No significant 
dependence on the parton distribution set was observed. 

A qualitative feature of the data, which contradicts the parton distributions studied, is the rapid decrease 
of d{x)/u{x) to unity beyond x = 0.2. At the same value of x as measured by NA51, x = 0.18, the d/u 
ratio is somewhat smaller than that obtained by NA51 Although the average value of Mf 



different for the two data sets, the effects of evolution on the ratio is expected to be small. 



IS 



To address the Gottfried sum rule violation, the extracted d/u ratio is used with the CTEQ4M value 



of d + u to obtain d — u at the average Mn+n- of the entire data set, M, 



7.35 GeV, 



25 



- d{x)/u{x) 
a[x) — u[x) = — 



d{x)/u{x) + 1 



u{x) + d{x) 



(55) 



The integral of d — u between x^in and 0.345 is calculated. Both the difference d — u and the integral 
are shown in Fig. ITO. The integral of the data is 



0.345 



dx[d{x) - u{x)] = 0.068 ± 0.007 ± 0.008 



(56) 



0.02 



compared to 0.076 and 0.1 for the CTEQ 4M and MRS (R2) sets integrated over the same region. When 
the range 10"'^ < x < 1 is considered, CTEQ 4M gives a value of 0.108 for the integral while the MRS 
(R2) resuh is 0.160. Above x = 0.345, both the CTEQ 4M and MRS (R2) distributions contribute 
less than 0.002 to the total integral. However, significant contributions to the integral from the parton 
distribution functions arise in the unmeasured region below x = 0.02, 30-50%, depending on the set. 
Such a large d/u asymmetry cannot arise from perturbative effects [^, especially since global fits to 
the parton distribution functions are tuned to accommodate the NMC, NA51, and, most recently, the 
E866 results {e.g. the MRST distributions [^] shown in the dotted curve in Fig. |]). 
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Fig. 10. Fermilab E866 values for (a) d — u and (b) the integral /^'^^^ {^d — tij dx' in the proton as a function 
of X. The curves represent the corresponding CTEQ 4M (solid) and MRS (R2) (dashed) results. The bar at 
0.147 lb 0.039 on the left axis in (b) shows the value obtained by NMC for the integral from to 1. Reproduced 



from Ref. |41]. Copyright 1998 by the American Physical Society. 



26 



A more complete picture of the individual contributions to — emerges from the NMC and E866 
data. Equation gives the decomposition of the difference into valence and sea components. Figure ^ 
shows the NMC data on the total difference as a function of x along with the E866 contribution to the 
difference, ocu — d. The curves show the total difference, the valence difference, oc u^ — dp, and the sea 
difference for two sets of parton distributions: MRS (R2), available before the E866 data, and MRST, 
which included the E866 data in the fit. Other relevant data in the MRST fit include an improved 
as determination from the CCFR data, the final NMC data on Fg*^, -F^'^, and F2^/F2^, and more 
precise information on the u/d ratio at large x from the Fermilab Tevatron [^. The effect of these 



improvements on the global analysis is clear from the plot. While both parton distributions fit the NMC 
data, only the MRST set agrees with the E866 data. Note that reducing the sea contribution to Eq. ( ^41) 
also significantly reduces the valence difference. 

The most recent experimental study of flavor asymmetries and the Gottfried sum rule is by the HERMES 
collaboration at HERA. They measured the charged pion yields in semi-inclusive DIS with a 27.5 
GeV positron beam on hydrogen, deuterium, and ^He gas jet targets in the kinematic range 0.02 < 
X < 0.3 and 1 < < 10 GeV^. Final-state hadron production requires the fragmentation function, 
D^{z), for parton i to produce hadron h with momentum fraction z of the initial parton momentum. 
The fragmentation functions Df [z) were extracted from the ^He data and assumed to be independent 
of the target. The HERMES collaboration determined the ratio 

d{x) - u{x) ^ J{z) [l-r{x,z)] - [l+r{x,z)] ^^^^ 
u{x) — d{x) J{z) [1 — r(a;, z)] + [1 + r{x, z)] 

where 

._ d(yf{x,z)/dz-dal{x,z)/dz 
^^^'^^ ~ da^^{x,z)/dz - da^'^{x,z)/dz 

'^'^-~5[l-DZ-i^)/Dfiz)) • ^''^ 

The ratio r{x, z) was determined from the hydrogen and deuterium data. No z dependence was observed 
in the measurement of the ratio in Eq. (0). Using the GRV 94 LO |^ parameterization of u{x) —d{x), 
they found 

0.3 

J c/x(S(x) = 0.107 ±0.021 ±0.017 , (60) 

0.02 

consistent with the E866 results, Eq. (|56|) , even though the E866 results are obtained at (M^+^-_) = 
(Q^) = 54 GeV^ while (Q^) = 2.3 GeV^ in HERMES, as shown in Fig. 0. The model calculations will 
be described later. 
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Fig. 11. The E866 difference d — u at {Q"^) = 54 GeV^ [^] is compared with predictions of meson cloud and 
chiral models. The HERMES data [^] are also shown. Reproduced from Ref. [42| with permission from the 



Annual Review of Nuclear and Particle Science, Volume 49 © 1999 by Annual Reviews. 

An experiment following up the E866 results but with a lower energy proton beam, 120 GeV, has been 



proposed at Fermilab |5^. This experiment would be able to increase the statistics at x > 0.1 and 
obtain results at larger x than available with the higher energy beam. 



Models which might account for SU(2) flavor asymmetry include meson clouds, Pauli blocking, chiral 
models, and charge symmetry breaking. These models will be discussed in the following sections. For a 



more extensive review of Gottfried sum rule violating mechanisms, see Ref. [56 . 



3.1.2 Meson cloud models 

Clouds of virtual mesons were shown to play a role in the dynamics of chiral models of nucleon structure 
such as the cloudy bag. In these models, the proton wavefunction can fluctuate into e.g. vr+ri and n^p. 
These virtual pion states have been used to explain the A decay width and the charge square radius of 



the neutron [RO]. While the proton wavefunction is in the virtual ttN state, the valence quarks and anti- 



quarks of the virtual state can contribute to the proton parton distributions. The original model |£7 
has been developed further recently to include a number of meson-baryon fluctuations with different 
form factors and meson-baryon-nucleon couplings at each vertex. Applications not only include the 
d/u asymmetry but also quark and anti-quark asymmetries between strange and heavier quarks. We 
introduce the model here and will discuss other applications in later sections. 
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The amplitude of the SuUivan process is obtained from the DIS amphtude by replacing the 7*p vertex 
in the DIS formalism by e.g. a nNN vertex 



M^=<X\eJ^{0)\n > — 



1 



(61) 



including the form factor, F^attvI^), coupling constant, ^f^^^, and n — N isospin coupling (p*-T. Replacing 
p by and M with in Eq. (|^), one obtains 



IQ^2 I yy 2 




-tdt 



{t - ml) 



nNN 



X 



:-f\' dyF^ - fMB{y) 



(62) 



where = —m^y^ /{I — y). Thus F2'^ is the convolution of the pion structure function with the light- 
cone momentum distribution of the pion, fuBiy) ll58| , |59| , |60| , |6Tl , |62| , normahzed by the isospin coupling. 
The same formalism is used for quark and anti-quark distributions. 

The direct interaction of a photon with the meson cloud surrounding a nucleon does not contribute to 
the Gottfried sum. However, it does not necessarily follow that the mesons do not contribute to the 
Gottfried sum and the u — d distribution, as discussed below. Two general approaches have been used 

30 ,^,62 1 the interaction of the recoil baryon with the virtual 



depending on whether ||58| , p9| or not 
photon is included. Both are in principle the same. The main differences are in the number of meson 
and baryon states included. 



The largest number of contributions has been considered in Ref. |^9|. They include ttN , pN, uN, aN, 
rjN, ttA, pA, KA, K*A, KT,, K*T,, KY*, and K*Y* states couphng to the nucleon. In most models 
however, only the ttN and vrA |^| nucleon couplings are included. Since the pions alone do not 
contribute to the Gottfried sum rule directly, the nNN and ttA^A contributions to the sum rule can be 



written as 58 



F^ix) - F^ix) 



l-x 



X , 



Z\-[u^{x)-d^{x)] 



(63) 
(64) 



where y' = x/{l — y) and fnN(A){y) are meson-baryon light-cone momentum distributions. The valence 
normalization Z is defined as Z = 1/(1 -|- N.^^ + A^r) where the probability of finding a pion in e.g. the 
ttA^A^ state is A^^ = /q dyf^N{y)- 
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The light-cone momentum distribution of the virtual meson in the nucleon is the sum over the meson- 
baryon momentum distributions 



fM{y) = T.fMB{y) , (65) 

B 

with 

f I \ 9mNB f -i.-L\t','rnN.,mB) rj-, , ..2 /aa\ 
fMBiy) = ^yjdt _ p [FMNBit]] (66) 

— 00 

as defined from Eq. (|62|). The meson-baryon momentum distributions fMB{y) are the probabilities of 
finding meson M in an MB configuration in the nucleon with fraction y of the nucleon momentum in 
the infinite momentum frame. These distributions must be related to the baryon distributions in the 
nucleon in the same configuration, fsAiiy), by 

fMBiy) = fBMii-y) (67) 

to conserve probability |^. The t-dependent function I{t,mN,mB) has two forms, 

t + {mB-mNy for 5 G 8 



12 m 'j^ m 'g 



for 5 G 10 



(6J 



depending whether the baryon B is part of the octet or the decuplet. The upper limit of the integral 
when B ^ N is given by 

2 

tm = mNy- . (69) 

Several functional forms, including monopole, dipole, and exponential, have been assumed for the meson- 
baryon vertex form factor. 



MNB 



monopole 



a 



dipole (70) 

(t-m|,)/A2 



= exponential . 



It has recently been pointed out that in the covariant formulation of Eq. (|66D with t-dependent form 
factors, the probability conservation of Eq. ( p7D could not always be achieved However, if the 
calculation is done with light-cone wavef unctions, e.g. the dipole form of FMNBif) above is replaced by 
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[(A^ + m|^)/(A^ + Sttat)]^ where s^at is the square of the pion-nucleon center of mass energy, probabihty 
conservation is satisfied |^ . 



Studies have been made of the hardness of the form factor. A hard form factor with a monopole 
cutoff 1.0 < Am < 1-4 GeV is needed to explain the deuteron D-state admixture and nucleon-nucleon 



scattering [56| although more recent A^A^ potential models allow for softer cutoffs, A^ ~ 0.8 GeV 



Softer cutoffs are also generally obtained in quark models, A^ ~ 0.6 GeV in the cloudy-bag model 1^5 
and 0.7 < A^ < 1-0 GeV in flux-tube type models |]66[. An even softer cutoff, A^, < 0.5 GeV p7 



arises 



from analyses of the flavor asymmetric distribution {u + d)/2 — s. However, later analyses of this flavor 
asymmetric distribution with more recent data resulted in a slightly larger cutoff, A^ 0.6 GeV [BT 



which could be consistent with the quark model estimates. As more mesons and baryons are added to 
the calculation, the cutoff becomes larger, A^ ~ 0.74 GeV, because the probability to find the bare 
nucleon is reduced due to the presence of the meson clouds. The dipole and exponential cutoffs are 
related to the monopole cutoff by A^ = 0.62Ad = 0.78Ae |l6T |. 



Mesonic contributions to an anti-quark distribution in the nucleon are given by the convolution of the 
corresponding valence anti-quark distribution in a meson with the light-cone momentum distribution 
of the meson in the MB state. The contributions are given by the equation 



1 



X q^ix, Q^) = J2 "MB / dy fMsiy) - Qm ( 1 , (71) 

MB I y \y J 

where the summations are over combinations of meson M = {it, K) and baryon B = {N, A, A, S, S*) 
states and apfB are spin-flavor SU(6) Clebsch-Gordon coefficients. Equation ( [7I| ) corresponds to Eq. 



(11) for 



From Eq. (|5^), the Gottfried sum becomes Iq = {Z/3){1 — N^^/S + 5A^/3). According to this equation, 
the failure of the sum rule is due not to the photon interaction with the virtual pion but to the interaction 
with the recoil baryons. 

When contributions from light meson and baryon states beyond vrA^ and vrA are included in Eq. (0), 
the Gottfried integral becomes 

1 } _ _ 

{Z + Ai) , with Ai= dx {ui + Ui-di- di)sn\\ (72) 



3 







where are the x integrated meson and baryon contributions to the quark and anti-quark distributions 
in the nucleons. These results, given in Table 0, should be multiplied by the probabilities of finding 
the meson-baryon states in the nucleon. Although the nucleon "core" satisfies the valence quark sum 
/ dx{u^ — d^) = 1/3, its probability is reduced by the normalization factor Z due to the presence of the 
virtual MB states. In this approach, no contribution from pions or rho mesons enters the sum and the 
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violation comes from the normalization factor Z and the baryon contributions because u + u — d — d 
vanishes in the pion. 
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Table 3 

Coefficients Ai and Bi in two different descriptions Reproduced from Ref. |5f:] with permission from 
Elsevier Science. 

The sum could be written in a different form to satisfy the valence quark sum exactly. Part of the 
meson and baryon contributions can be factored into the valence sum, 1/3, and the deviation from 1/3 
is identified with the flavor asymmetry due to the Sullivan processes from the meson cloud, / dx{u—d)^^^\\ 
The Gottfried integral then becomes 

1 } _ _ 

Ig = -{l + Y^Bi) , with Bi = J dx {ui - rfi)suii • (73) 

* 

The coefficients Bi in Eq. ( |73D are also given in Table ^ The virtual ttB states contribute to the 
renormalization of the valence-quark distributions. Therefore the pionic renormalization contributions 
may be factored into the valence quark sum factor in Eq. (|73|) |^. Then the pion contributes to the 
deviation from the Gottfried sum. Because flavor symmetry is assumed in the MB states, only the pions 
and rho mesons contribute to the violation. Thus the two different mesonic descriptions are equivalent. 

Although there are slight model differences, the meson-cloud approach is successful in explaining the 
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major part of the NMC and E866 results. The sohd curve in Fig. |TI] shows a calculation with the 
meson-cloud model with ttNN and ttNA contributions. A dipole form factor is used with = 1 GeV 
for the irNN coupling and 0.8 GeV for the ttNA coupling, corresponding to A^ = 0.62 GeV and 0.5 



GeV respectively. Reasonably good agreement with the E866 d — u data is found |^ 



Even though meson cloud models, such as the calculation shown in Fig. |TT], may successfully explain the 
difference d — u, they usually have trouble reproducing the ratio d/u. An approach that has had relative 
success in obtaining agreement with both the difference and the ratio is the inclusion of a perturbative 

based on the assumption that the proton 



component in the meson cloud p8|,|69[]. The analysis of Ref. pS 



is composed of three dressed quarks or valons in the infinite momentum frame |70|, is discussed in more 



detail here. One of the valons may emit a gluon which subsequently decays into a qq pair. This qq pair 
can then nonperturbatively rearrange itself with the valons of the proton into a meson-baryon bound 



state. The flavor symmetric GRV distributions |^ are assumed to be the underlying perturbative sea. 



The ratio d/u is computed assuming two components in the sea, the nonperturbative meson cloud and 
a perturbative contribution, 



-;NP 



djx.Q^) ^ d ix,Q^)+tix,Q^ 



(74) 



The difference d — u is extracted from the d/u ratio according to Eq. (pq) , as also done in the E866 
analysis [^, rather than from the difference between the numerator and denominator of Eq. (fTil). The 
perturbative anti-quark distribution in Eq. ([7^) is taken from the GRV LO symmetric sea |]7l| . 



The nonperturbative u and d distributions are calculated in the meson cloud model assuming that the 
d is generated from vr+n and 7r+A° fluctuations while the u arises from tt^A^~^ fluctuations. 



-NP / ^2 N /' 



y L 








(75) 
(76) 



The fluctuations containing 7r° mesons, e.g. tt'^A^, are not included because flavor singlets such as 
dd — uu are assumed to be more likely to annihilate than to recombine with the valons in the 



proton. The c?^ and 11^^ valence distributions are taken from the GRV LO pion distributions |7^. The 
pion probability densities from ttB states are calculated by the convolution of a recombination function 
associated with pion formation, R{x, y, z), 



Rix,y,z) 



7 — d[x 

X 



y 



(77) 
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and a valon quark distribution function, F{y^ z), 



F{x, y) = (3yv{y)zq{z){l - y - z)" 
so that 



PM^) =1—1 —F{y,z)R{x,y,z) . (79) 
J y J z 

In Eq. (0), the exponent a is fixed by requiring that the pion and baryon in the fiuctuation have the 
same velocity, resulting in a = 12.9 for 7T~^n and 18 for the ttA fiuctuations. The quark distribution is 
calculated from the decays of the valon to gluons, v ^ v + g and gluon splitting to qq, g q + q, 

= Ql) = j '^p„ I (f) (80) 

X y 

where Qy = 0.64 GeV^ and P^g and Pgg are the valon and gluon splitting functions for QCD 

evolution, as in Eqs. (p^) and (|37|). The valon distribution, independent of Q^, is v{x) = (105/16)-y^(l — 
xy. The normalization constants N, 7, and /3 are fixed from the difference d — u in Eq. ([551) . 



The results for the ratio, Eq. ([7^), and the difference, Eq. (pBD, are shown in Fig. The solid curve is 
the result at the valon scale = 0.64 GeV^ while the dashed line is a pseudo-evolved result at = 54 
GeV^ Both calculations are in reasonable agreement with the measurements due to the presence 



of the perturbative component. The calculation describes the ratio as well as the difference because the 
valon model fixes the normalization of the perturbative component while the meson cloud model shifts 
the contributions to d and u between the meson and baryon components. Thus the meson cloud model 
alone can only describe the difference. Even better agreement is obtained if the valence distribution in 
the pion in Eqs. ( |75D and ( [7B| ) is multiplied by an additional power of x [^]. This modified pion valence 
distribution is not excluded because the low x pion parton distributions are not well known. 



3.1.3 Effects of Pauli blocking 

The Pauli exclusion principle forbids two fermions from having identical quantum numbers. Thus, if 
a quark in the nucleon wavefunction has a certain set of quantum numbers, a g in a virtual qq state 
cannot be generated in the sea if it has the same quantum numbers as an already existing q. The new q 
is Pauli blocked. Because the proton has two valence u quarks and one valence d quark, uu pair creation 
is more likely to be Pauli-blocked than the dd pair, resulting in a difference between u and d, as first 
suggested by Field and Feynman ||7^ . 

Pauli blocking effects in the nucleon have been calculated in 1 + 1 dimensions using the MIT bag model 
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Fig. 12. The E866 ratio d/u and difference d — u are compared to calculations at = 54 GeV^ |41] with the 
meson cloud model of Ref. [^]. The solid curves show the and d distributions without evolution 
(solid curves) and with pseudo-evolution (dashed curves). Reproduced from Ref. [38|. Copyright 2000 by the 
American Physical Society. 



58| , [75| . In 1 + 1 dimensions, there are three color states for each flavor. In the ground state, two of the 
three w-quark states and one of the three d-quark states are occupied. Only one more u quark can be 
added to the ground state, but two more d quarks can be accommodated. Since the additional u or 
d quarks in the ground state would be accompanied by u and d anti-quarks, the expected sea quark 
asymmetry is rather large, d = 2u. 



No 3 + 1 dimensional calculation exists at this time so that the Pauli blocking effects are obtained from 
a naive counting estimate. In the 3 + 1 dimensional case, there are six states, three color and two spin, 
in the ground state. There are four available ground states for u quarks and five states for d quarks, 
so that d = 5u/4:. Because no valence anti-quarks exist in the bag, the sea quark contribution comes 
from a quark being inserted, interacting in the bag, and then being removed. Therefore, the d excess is 
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related to the distribution of a four-quark intermediate state f^i^x) 



dx [d{x) — u{x)\ = / dxf4{x 



where P2 is the integral of a distribution associated with a two-quark intermediate state. Because the 
u and d distributions have not been calculated, the Pauli contributions were assumed to be x^^fl — x)'^ 



in Ref. 

used, and a large exponent 



The constant A is chosen to match the small x behavior of the valence distributions 
B = 7, was chosen so that Pauli blocking contributes only at small x. 

The overall normalization is not 



This value of B is consistent with counting rules for sea quarks 
determined theoretically but 1 — P2 is expected to be ~ 0.25 because of the naive counting estimate in 
3-1-1 dimensions. 



Pauli blocking effects could in principal, produce the d over u excess. Unfortunately, the qualitative x 
dependence of the effect is only calculated in 1 + 1 dimensions. It was suggested that a 10% Pauli-blo eking 
effect combined with pionic contributions from the meson cloud could explain the difference — F2"' 



58| and the E866 d/u ratio and d — u difference relatively well [^. However, if antisymmetrization of 
the quark wavefunctions is included, the conclusion could change dramatically |^ . The same m- valence 
excess which suppresses the uu pair creation produces extra uu pairs in the antisymmetrization of the 
second u. This additional source of u could cause m > c?, in contradiction to the NMC and E866 results 



3.1.4 Chiral models 



Several different types of chiral models have been proposed to explain the Gottfried sum rule violation. In 
a chiral field theory with quarks, gluons, and Goldstone bosons (pions) [^,78,7S,BC], flavor asymmetries 



arise from virtual photon interactions with the pions. On the other hand, in chiral soliton models 



3 1,32, 83 1 a portion of the isospin of the nucleon is carried by pions. The Gottfried sum rule violation 



is then determined from the moments of inertia of the nucleons and pions. An extension of the linear 
sigma model to the Gottfried sum rule has also been proposed p4|,|85|,p6| . In this section, only the chiral 



field theory approach is briefly described. For more details, see Ref. |56 



Because chiral symmetry is spontaneously broken, any description of low energy hadron properties 



should include this symmetry breaking. The effective Lagrangian in chiral field theory is |56 



f82) 



where ip is the quark field and is the covariant derivative. The vector and axial-vector currents are 
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defined in terms of Goldstone bosons so tliat 



(83) 



wliere ^ = e'^'^^/^'^ and tlie +-sign goes witli tlie vector current. Tlie meson matrix 11 = g (8> g is 
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7r7V2 + We 



TT 



K- 



\ 



\ 



TT 

K- 
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S4) 



Expanding V^^ and in powers of H// gives = 0(11//)^ and = id^Ii/ f + 0(11//)^. Then the 
quark-boson interaction, proportional to A^, becomes £ng = —{.9A/f)i^^tJi^'l^l^''P- 

The Gottfried sum rule violation arises from q —>■ splittings where the final-state pion is massless. 
The u quark can split into either u — *• 7r~^d —>■ udd or m — >• tt^m — u{uu), u{dd). Thus if a is the splitting 
probability for u vr+c?, 



M — s> avr -|- ad H — vr H — -u = -- Tm + 5a + -u + 5a 



^5) 



Likewise, 



_ a Q a , a r 

d an + au -\ — tt H — a 



5m + 7(i + 5m + d 



2 2 4 

The probability a is determined from the splitting function, Png'^q, defined as 



^6) 



i-n.'^,(^) - ^^^^^ z 



dt 



A2 



2N2 



(87) 



where, as in the meson-cloud model, = m'^z — ■m'^,z/{l — z), similar to Eq. (|5^), and A :^1.169 GeV. 
This large value of A is assumed to be the scale at which chiral symmetry is broken. Integrating over t 
for the process u 7c~^d gives 



gl ml 



J dze{A^ -T{z))zl\n 
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h? + Ml t{z) + Ml 
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where 9{x) = 1 when x > and zero otherwise. Since ^ m^, t(z) = 
With A ^ 1.169 GeV, the Gottfried integral becomes = {I - 2a)/3 
Ig- 



0.278. Increasing A decreases 



The anti-quark distribution as a function of x is |177| , [78 



j,k,i \ 



6ihSii\'^ } dy j dz _(u) f X 



Jjkfil 



y J z 

X x/y 



Puk^j{z)q^j{y) 



^9) 



where the indices j, k, and / are summed over flavor. The d — u difference d — u is calculated at Q'^ = 54 
GeV^ using the formulation of Szczurek et al. [^| and compared with the E866 and HERMES data in 
Fig. 0. The difference in the chiral approach is concentrated in the small x region, even more so than 
in the meson-cloud model, as seen in the dashed curve of Fig. 0. This concentration at low x occurs 
because pions are coupled to constituent quarks in the chiral model and the constituent quarks carry a 
smaller fraction of the nucleon momentum than the virtual pions in the meson cloud. 



3.1.5 Charge Symmetry Violation 



Charge symmetry, Eq. (j2^), is respected to a high degree of precision in nuclear physics with most 
low energy tests finding that it holds to ~ 1% in reaction amplitudes [^. Thus, charge symmetry is 



generally assumed to be valid in strong interactions and all global analyses of parton distributions [ P5| . 

Experimental verification of charge symmetry is difficult for two reasons: any violation effects are ex- 
pected to be small and charge symmetry violation may mix with or be misinterpreted as flavor symmetry 
violation. Thus, Ma suggested that the interpretation of the Gottfried sum rule violation and the 



Drell-Yan ratio as flavor symmetry violation could instead be the result of charge symmetry viola- 
tion. Given such ambiguities, experiments that could distinguish between charge and flavor symmetry 
violation are needed. Some have already been proposed 



Possible effects of charge symmetry violation have recently been examined by Boros et al. They 
define charge symmetry violating ,CSV, distributions 

5u{x) = Up{x) — dn{x) 6d{x) = dp{x) — Un{x) (90) 
6u{x) = Up{x) — dn{x) 6d{x) = dp{x) — m„(x) , (91) 

which would disappear if charge symmetry holds. The isoscalar structure functions defined in Eqs. (I 



(|33|) are then modified to include the violation pT 



F2csv{x) = - Tsl'^iHx) + ^d{x)) + 5u{x) + 5u{x)] (92) 
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(96) 



Likewise, one can form a charge ratio relating the neutrino structure function to the charged lepton 
structure function similar to Eq. (|39|) 
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(97) 



where Q(x) = J2q=u,d,s{lix) + q{x)) — |(s(a;) + s{x)). The denominator of Eq. (p7|) was expanded 
retaining only leading terms in the small quantities s{x) — s{x) and Sq{x). If Rc{x) 7^ 1 at any x, then 
either s{x) 7^ s{x) or charge symmetry is violated. The CCFR collaboration compared their neutrino 
structure function F^^^"^ |22] with measured by the NMC collaboration [13| and found that for 
x < 0.1, i?c ~ 0.9 in the region of overlap between the two experiments. When nuclear shadowing 
corrections were applied to the targets, Rc remained inconsistent with unity. It therefore appears that 
in this region, s{x) 7^ s{x) or charge symmetry is violated. 

The CCFR data allow s{x) and s{x) to be determined independently. Such an analysis by the CCFR 



collaboration ||2^ will be discussed in more detail in Section 3.2.1. Here, the CCFR dimuon data 
were used to set limits on charge symmetry violation effects [^]. Dimuon production from i/^ and 77^ 
interactions with the strange nucleon sea is defined as 



xs 



{x) = -x [s{x) + l{x)] + -{213' - 1) X [s{x) - s(x)] 



(9J 



where experimentally (3' = NI^''/{NI^'' + N^^) ^ /5 in Eq. (g|). Combining Eq. (|8]) with Eq. (|9|), one 
can simultaneously solve for s{x) and s(x). The anti-strange distribution determined in this fashion is 
negative and thus unphysical [B^. If charge symmetry is assumed to be violated, the difference between 



Eqs. (|3|) and 
5 



6 



''yx)-?>Ff\x 



X 



) can be used to estimate the size of the violation [ pT[] 

x{2(3-l 



xs^^ix) 



[s{x) - s(x)] 



+ - [(5/5 - l){5d{x) - 5u{x)) + (4 - hP){5d{x) - 6u{x))] 
6 



x(2/3 - 1) 



s{x) — s{x)] + 2 ^ [5d{x) — 5u{x)] . 



(99) 



39 



In the last line, it is assumed that (3 = (3' and 6q^ = 6q — 5q ^ 0. The assumption 6q^ ~ means that 
charge symmetry violation should manifest itself most strongly in the sea. Fig. show the magnitude 



of the effect obtained for the two extremes s{x) = s{x) and s{x) = [^. The result is surprisingly 
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Fig. 13. Degree of charge symmetry violation obtained from Eq. (^) assuming that s{x) = s{x) (open circles) 
and s{x) = (filled circles). Statistical and systematic errors have been added in quadrature. Reproduced from 
Ref. pH]. Copyright 1998 by the American Physical Society. 



6u{x) [91] which has the effect that Up + dp = Un + dn- The results of 



large. Theoretically, 5d{x) ^ 
Fig. ^ suggest a charge symmetry violation of ~ 25% at low x. Such large effects, if they exist, would 
necessarily have to be incorporated into future fits of parton distribution functions. 

This rather large charge symmetry violation at low x would be mixed with flavor asymmetry in the DIS 
and Drell-Yan studies of the Gottfried sum rule violation. In this case, with —5u ~ 5d^ Eq. (0) would 
be rewritten as 



X 



IP I 



X] 



rrix)) 



\ K(x) - d'^ix)] + I [u{x) - d{x) + 5d] . 



(100) 



Since measurements of Iq suggest that the second term on the right hand side of Eq. (|100|) is negative and 
5d is positive in Fig. |TB|, charge symmetry violation would still require a large SU(2) flavor asymmetry, 
even larger than that suggested without charge symmetry violation. Note that the magnitude of Iq 
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would not change, only the interpretation of the result if charge symmetry violation mixes with flavor 
asymmetry The d/u ratio with and without charge symmetry violating effects is shown in Fig. |n. 
The open points show how the flavor asymmetry would need to be enhanced for Iq to remain the same 
after charge symmetry violation is accounted for. This would indeed require a substantial modification 
of the low X d and u distributions in the global analyses of parton distribution functions. 




Fig. 14. The ratio d/u as a function of x for the NMC [^] (triangles) and E866 [^T| (circles) data. The solid 
points do not assume any charge symmetry violation while the open points include charge symmetry violation 



at low X as parameterized in Fig. 13. Reproduced from Ref. [9^]. Copyright 1999 by the American Physical 
Society. 



Another test of charge symmetry violation can be made with production at hadron colliders. Boros 
et al. suggest comparing and W~ production in pD collisions at the Relativistic Heavy Ion 
Collider ||93| at Brookhaven National Laboratory, BNL, and the Large Hadron Collider at CERN. 
However, as recently pointed out by Bodek et al. the DIS and Drell-Yan measurements discussed 
here cannot extract charge symmetry violating effects alone. Measurements of dp/up or dn/un would 
isolate these effects more definitively. At the Fermilab Tevatron, and W~ are produced by ud 
and du — > W~ respectively. The u quarks carry more of the proton momentum than the d quarks, see 
Figs. so that the follows the incoming proton while the W~ follows the anti-proton, creating a 
charge asymmetry between the two in rapidity. This asymmetry, measured through the decay leptons. 
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is directly sensitive to the dp/up ratio Bodek et al. parameterized the charge symmetry violation 
of Fig. in two different ways. In the first, the average of the d and u distributions is unchanged by 
charge symmetry violation while in the second, the entire violation is placed on the d, leading to a shift 
in the total sea quark distribution. Both parameterizations are compared to the Tevatron CDF data 
at 1.8 TeV in Fig. |15|. The large charge symmetry violation suggested in Refs. [^,0 overpredicts 
the asymmetry for y > while calculations with parton distribution functions assuming no charge 
symmetry violation agree with the data. Bodek et al. proposed that even though there is no evidence 
for charge symmetry violation in the CDF data, the large violation proposed by Boros et al. pT| , |92| can 



be preserved if charge symmetry is violated only in the neutron parton distributions. In this case, the 



neutron sea would be larger than the proton sea. Thus, while the results of Ref. weaken the case 
for charge symmetry violation, they cannot eliminate it until better measurements can be made. Note 



by Boros et al. [pT| , p2l is at leading order |9^. A treatment beyond 
leading order might eliminate the need for charge and flavor symmetry violations in the neutron. 



however that the treatment of 



0.3 



6 

S 0.2 

>^ 

< 

bp 0.1 







-0.1 



-0.2 







Model 2 



Model 1 




CTEQ4M (CSV = 0) f 



CTEQ4M (CSV = 0,MOD) 



CDF 



0.5 



1.5 2 

Lepton Rapidity 



Fig. 15. The CDF asymmetry as a function of rapidity [^ ] compared to calculations without charge 
symmetry violation, labeled CTEQ4M (CSV=0) and CTEQ4M (CSV=0,MOD), and with the effect of Fig. |T3 
using the two parameterizations of Bodek et al. |^, labeled Model 1 and Model 2. Reproduced from Ref. p5| . 
Copyright 1999 by the American Physical Society. 
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3.1.6 Light sea summary 



Experiments have clearly shown that not only isu^d but also that the difference is a strong function 
of X p9| , |il|] . A number of models have been proposed to describe this difference, including meson clouds, 
Pauli blocking, chiral models, and charge symmetry violation. All but the last assume that u ^ d \s 
due to SU(2) flavor asymmetry. Meson cloud models have generally had the most success describing 
the difference d — u but, with the exception of Ref. generally have not had as much success with 
the ratio d/u. Chiral models tend to concentrate the difference at low x, underestimating d — u as a. 
function of x. The proposed charge symmetry violation would also require a large flavor asymmetry 



21,921 but relies on the accuracy of the CCFR strange quark measurements 22 



3.2 The strange sea 



Possible manifestations of nonperturbative effects in the strange sea are now discussed along with 
some phenomenological explanations. The CCFR data which might imply a difference between the s 
and s distributions ||2^ are discussed first. Three experiments which indirectly measure strange form 



factors, two low-energy elastic ep scattering determinations of the parity-violating asymmetry [|97| , |98 
and neutral current v and a V scattering experiment on nucleons 



99 , are also introduced. 



3. 2. 1 Experiments 

The CCFR collaboration has performed LO and NLO evaluations of dilepton production in neutrino 
scattering off nucleons to more directly study the strange quark distributions [^J2^. Although, as 
mentioned earlier, they found s{x) = s{x) within their uncertainties, the result is not inconsistent with 
s{x) 7^ s{x) in some x regions, as is now discussed. 

Pairs of opposite sign muons can be produced via charm production in and scattering off nucleons: 
i/^ + AT — , + X (101) 

+ N — ^/i++c + X (102) 



In the first case, the neutrino interacts with an s or d quark, iJ + s,d c+j2~, to produce a charm quark 
which then decays semileptonically to a Because the ci-induced channel is Cabibbo suppressed rela- 
tive to that of the strange quark, this process can provide a determination of the strange sea. Likewise, 
the anti-neutrino-induced interaction is an independent measurement of the anti-strange distribution. 
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The non-strange quark and anti-quark components of the sea were assumed to be symmetric so that 
u{x) = Us{x) and d{x) = ds{x). An isoscalar correction was also apphed assuming u{x) = d{x). The 
strange sea was allowed to have a different magnitude and shape than the non-strange sea. Shape 
parameters were determined from two fits, both starting from the same set of proton parton distributions 
at initial scale /Xq = 1 GeV^ [[IOCI|| . In the first fit, s{x) = s{x), while in the second s{x) and s{x) were 



determined independently. The x-integrated strange fraction of the sea is set by the parameter 



K 



Jq dx [xs{x, /i^) -|- xs{x, yU^)] 
Jq dx [xu{x, /i^) -|- xd{x, /i^)] 



(103) 



where k = 1 indicates an SU(3) flavor symmetric sea. The shape of the strange quark distribution is 
related to that of the non-strange sea by a shape parameter S. li S = 0, the strange sea would have the 
same x dependence as the non-strange sea. 



Assuming that s{x) = s{x), the sea quark distributions were parameterized as 

xu{x, ji^) + xd{x^ jj^) 



xs{x,^ ) = As{l - x) 



xq{x,^i^) = {2 + As{l-xf) 



The normalization Ag in Eq. 
the parameters are 



xu{x, /i^) + xd{x, /i^) 



is defined in terms of k and 5. With the requirement s{x) 



- 0.051 - 0.017 
^ 0.050 + 0.036 
4- 0.66 + 0.08 

- 0.60 - 0.20 



K = 0.477 

5 = -0.02 
= 1.70 ± 0.19 ± 0.02 GeV . 



(104) 
(105) 
s(x), 

(106) 
(107) 
(108) 



Note that the charm quark mass enters the fit through the original charm production reaction in 
Eqs. (|101D and ([TO^) . 



The strange quark, Eq. (|10 
this fit are shown in Fig. IT 



and the total anti-quark, Eq. ( |105| ), xq = xu + xd + xs, distributions in 
at the initial scale /ig = 1 GeV^. The strange quark distribution is scaled up 



by a factor of 2.5 in Fig. |T6](a) to show its similarity to the anti-quark distribution. Since 6^0 would 
indicate a shape difference between xq{x) and xs{x), the fit value S = —0.02 1 ^go indicates no shape 
difference at NLO. In a previous leading order fit, they found that the strange quark distribution was 
softer than the total quark sea by (1 — a;)^'^^°'^. The difference between the LO and NLO fits may be 
attributable to a softer xq{x) at NLO than LO. This similarity is highlighted in Fig. |l6|(b) where the 
ratio of the strange quark distribution to xq is shown. The ratio is essentially constant over all x. 

If one allows the momentum distributions of the s and s quarks to be different, the sea quark distribu- 
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Fig. 16. The CCFR xq and xs distributions in the s{x) = s{x) fit are compared at the initial scale 1 GeV^. 
In (a) the individual distributions arc shown with the strange quark distribution (solid curve) multiplied by a 
factor of 2.5. In (b), the ratio xs/xq is shown. 

tions can be parameterized by: 



xs{x, ji^) = As{l — x) 



xs{x, jj^) = A'g{l — x) 



xu{x, 11^) + xd{x, iJ?) 



S-AS 



xu{x, jjj^) + xd{x, /i^) 



xq{x,i^'') = {2 + -{A,{l-xY 



xu{x, /i^) + xd{x, /x^) 



(109) 
(110) 

(111) 



The fit required the number of s and s quarks to be equal: 



1 ± 
J dx s{x, fjL^) — J dxs{x,^^) 



(112) 





The normalization factors Ag and A'^ are defined in terms of 5 and A5. In this case they found 
K = 0.536 ± 0.030 ± 0.036 ; ± 0.009 



5 = -0.78 ± 0.40 ± 0.56 ± 0.98 ± 0.50 
= -0.46 ± 0.42 ± 0.36 ± 0.65 ± 0.17 
rur = 1.66 ± 0.16 ± 0.07 t ± 0.01 GeV . 



(113) 
(114) 
(115) 
(116) 



In this analysis, the total strange quark content, given by is slightly higher and the shape of the 
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s quark distribution is different from the non-strange sea. Compare Eqs. ( |113| ) and (|114|) with the 
values of k and 6 in Eqs. ( |10(j| ) and ( |107| ). The charm quark mass parameters in the two fits are very 
similar. The uncertainty in A6 is large enough for the s and s momentum distributions to be consistent 
with a difference limited to —1.9 < A6 < 1.0 at the 90% confidence level. The resulting s{x) and 
s{x) distributions in Eqs. ( |109| ) and ( |110| ) are compared to the q distribution, Eq. ( |111| ), in Fig. |17 



at the initial scale /Iq = 1 GeV^. Both the s and s distributions are multiplied by a factor of 2.5 
to facilitate comparison with the larger q distribution. The strange quark distribution is larger than 
the anti-strange distribution and deviates more from the total anti-quark distribution than does the 
anti-strange distribution. 
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Fig. 17. The CCFR xq, xs, and x's distributions in the s{x) ^ s(x) fit are compared at the initial scale 1 
GeV^. In (a) the individual distributions are shown with the strange quark (solid curve) and anti-strange 
quark (dot-dashed curve) multiplied by a factor of 2.5. In (b), the ratios xs/xq (solid curve) and xs/xq (dashed 
curve) are shown. 



It is important to note that the starting distributions ||100|| used in the fits assume fiat sea and gluon 
distributions a.s x —>■ and a flavor symmetric light quark sea. Both these assumptions have had to be 
modified recently. If the analysis was repeated with newer starting distributions which take these new 
developments into account, perhaps the results would change. In particular, the discrepancy between 
these data and the muon-induced DIS data which led Boros et al. to argue for a large charge symmetry 



violation |2T| could be reduced and the difference, if any, between s{x) and s{x) could be more clearly 
delineated. 



A step in this direction was recently taken by Barone et al. ||101|| in an analysis of DIS data that 
emphasized the u and V induced cross section data. The CCFR data were not used in their global fit 
but the fitted strange quark distribution was shown to agree with their data. Barone et al. also allowed 
s{x) 7^ s(x) in one of their fits to all DIS data but did not constrain the shapes to be proportional to 



0.5[m(x) -|-(i(x)] as in Eqs. ( |109| ) and ( |110|) . Since the DIS data included in the Barone et al. fit was more 
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balanced between v and v induced events with better statistics, tighter constraints on the difference 
between s{x) and 's{x) were set. They note that the assumption s[x) ^ gave a better overall fit to 



the data with xs{x) harder than xs{x). The difference is largest at x > 0.4, as seen in Fig. ^ flOl 





Fig. 18. Results of a fit to all DIS data allowing s{x) ^ [101|. The results for the difference xs{x) — xs{x) 
is shown in (a) while the ratio s{x)/s{x) of the fit is compared to the CCFR result at = 20 GeV^ in (b). 
Reproduced from Ref. 



101] with permission from Springer- Verlag. 



Strange form factors have been inferred from the parity violating asymmetry in electron-proton elastic 
scattering 

_ -GfQ' eGlGl + TGlfil, ' 0-5(1 - 4 sin^ d w)e'Gl,F^ ^^^^^ 



7ia 



V2 e{Gir + r{Giy 



where r = Q^/Am^, e = 1/(1 + 2(1 + r) tan2(6'/2)), e' = y'r(l + r)(l - e^), and 9 is the scattering 
angle. This asymmetry between right and left handed incident electrons is due to the interference of the 
electromagnetic and neutral weak amplitudes. The electromagnetic electric and magnetic form factors 



47 



are and G],/ while the neutral weak electric, magnetic, and axial form factors are Gf;, G^,j 
To leading order, the proton weak magnetic form factor is 



^ and Ff . 



4 sin^ e 



Wj 



G 



M 



118) 



The electroweak radiative corrections to Eq. (|118|) have been calculated ||102|] . The weak mixing angle 



6w has been determined with high precision and G^i and G]^ are known. 



The strange magnetic form factor has been determined by the SAMPLE collaboration at backward 
angles. A 200 MeV polarized electron beam is directed on a liquid hydrogen target. The polarized 
electron source is a GaAs photoemission crystal. The photoemission is stimulated by a laser beam 
incident on the crystal. The laser beam is circularly polarized and the electron beam helicity is reversed 
by changing the voltage to the cell which polarizes the laser beam, reversing the circular polarization 
of the light. The helicity is chosen at random for ten consecutive pulses and the complement helicities 
are used for the next ten pulses. The scattered electrons are detected by a large Cerenkov detector 
which accepts momentum transfers of ~ 0.1 GeV^. The asymmetry is found for pairs of pulses 
separated by 1/60 s to minimize systematic uncertainties. Each "pulse pair" asymmetry is equivalent 



to a measurement of the parity violating asymmetry, Eq. ( 117 ). 



The SAMPLE collaboration determined Gf^ from A. The strange magnetic form factor G\^ can be 
obtained from A since the backward scattering angles and low cause the second term in the numerator 
of Eq. ( |117| ) to dominate the asymmetry. They found 



A- 



: -6.34 ± 1.45 ±0.53 ppm , 
: 0.34 ± 0.09 ± 0.04 ± 0.05 n.m. 



(119) 
(120) 



at 
If 

corresponds to |9 



0.1 GeV^. The last uncertainty on Gfj is due to theoretical uncertainties in the form factors. 



If G%j = 0, then Gfj 



0.40 n.m. [^]. Thus the difference between the expected and the measured G 



M 



Gl^iQ^ = 0.1 GeV^) = 0.23 ±0.37 ±0.15 ±0.19 n.m. 



:i211 



consistent with zero. 



The HAPPEX collaboration performed a similar measurement of A at a higher energy, with 3.356 GeV 
electrons, and the scattered electrons were measured at more forward angles (6'iab) ~ ±12.3° [Q. The 
polarized source was a GaAs photocathode excited by a circularly polarized laser. The helicity was set 
every 33.3 ms and structured as pairs of consecutive 33.3 ms periods of opposite helicity. The presence of 
false asymmetry was ruled out by inserting a half-wave plate in the laser beam. The correlation between 
the half-wave plate and the sign of the asymmetry was an unambiguous signal of parity violation. The 
measured asymmetry was 
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A = -14.5±2.0±l.lppm (122) 

at = 0.48 GeV^. Since r and e are of the same order at forward angles, G% and G%i cannot be 
separated from each other. Thus the combination of the strange electric and magnetic form factors was 
reported [@ 

7 
M 



tG 



G% + -7^<^M = 0.023 ± 0.034 ± 0.022 ± 0mQ{5Gj) n.m. . (123) 
e Ge 

The last uncertainty is an estimate of the theoretical uncertainty on G^"". A new run should improve 
the precision of Eq. ( |123| ) by a factor of two ||103 |. 



Both the SAMPLE and HAPPEX results are consistent with strong suppression of the strange sea at 
low energies ||104|| . To separate G% from G%.^ at fixed Q^, a number of scattering angles and thus beam 



energies are needed |p.03|| . An isoscalar target would also help separate G% and G%j [^. Additional 



measurements of weak currents and strange form factors are planned at Mainz |105| and the Thomas 



Jefferson Nation Accelerator Facility, TJNAF [106 . 



Information about the strange form factors of the nucleon can also be obtained from scattering of v 
and V on nucleons. The BNL experiment E734 performed such an investigation using wideband 



neutrino and anti- neutrino beams on a ^^C target with an average energy of 1.3 GeV and momentum 
transfer 0.5 < < 1 GeV^. The majority, 80%, of the data were due to quasielastic proton knockout 
from the carbon nuclei, vn fi'p and Vp fJ'~^n, charged current processes, and the remaining 20% 
arose from u and 17 elastic scattering on free protons, up up and up up, governed by neutral 
current interactions. Since the quasielastic interactions do not involve strange quarks, the strange form 
factors only appear in neutral current processes. 

Experimental differential cross sections are obtained by a convolution of the expected differential dis- 
tribution with the neutrino/anti- neutrino energy spectrum, (f)„{eu), so that 

, , \ f , da,, , ^ , ^ 



■=111111 



where emin = 0.2 GeV and emax = 5 GeV. In Eq. ( p.24| ), u represents either neutrinos or anti-neutrinos. 
The total neutrino flux is In the elastic scattering events, is directly obtained from the kinetic 
energy of the scattered proton, Tp, in the laboratory frame, = 2mpTp. In the quasielastic knockout 
events, an effective momentum transfer for u and u scattering off free nucleons was obtained by correcting 
for the binding energy and Fermi motion of the struck nucleon. The total cross sections, {ciy), were 
obtained by integrating Eq. ( |124|) over in the range of the data. The cross section ratios were |]99| : 



49 



R = ^^"P^'^p'^ = 0.302 ± 0.019 ± 0.037 (125) 

i?^ = -^^^^^=^ = 0.153 ±0.007 ±0.017 (126) 

i?57=y^^^^^ = 0.218 ±0.012 ±0.023 (127) 

where R is the elastic u/u scattering ratio while Rjy and R^j are the ratios of elastic scattering to 
quasielastic knockout cross sections after the correction to "free scattering" . Thus R^ and Rtj are ratios 
of charged to neutral current processes while R involves only neutral current scattering. 

In the following sections, models of the interactions discussed here are introduced and compared to the 
data. 



3.2.2 s/s asymmetry 



The CCFR NLO evaluation of the strange quark distributions did not exclude the possibility that the 
strange and anti-strange quark distributions are different. The analysis indicated that this asymmetry 
is small, on the order of (1 — x)^°-^^ [Q, as suggested by Eqs. ( |109|) , (|110|) , and ( |115| ). There is actually 
no reason why the sea quark and anti-quark distributions must be the same besides phenomenological 
prejudices as long as their total numbers are identical. Nothing in the QCD Lagrangian explicitly relates 
the sea quark and anti-quark distributions in the nucleon. Charge conjugation symmetry only tells us 
that the quark distribution in a proton is the same as the anti-quark distribution in an anti-proton 



El 



Ji and Tang showed that for every sea quark interaction, the corresponding sea anti-quark interaction 
may be obtained by changing the direction of quark propagation in the Feynman diagram ||107|| . However, 
different color factors are associated with the quark line direction, changing the relative interaction 
strengths. They found that while the quark and anti-quark diagrams are of the same order in l/N^ {N^ 
is the number of colors) the numerical coefficients are not constrained to be identical. 

A difference between the strange and anti-strange quark distributions can arise naturally in the meson 
cloud-type models introduced in Section 3.1.2. Since the proton can fluctuate into a virtual kaon-hyperon 
intermediate state, the strange anti-quark distribution in the kaon is different from the strange quark 
distribution in the hyperon ||108|| . The results of meson cloud model calculations [|109|J110| , |111| , |112|| will 
be discussed in this section followed by a calculation of the strange and anti-strange distributions based 
on the locality of the strange sea ||107| . 



The first model, by Paiva et al. |109|, follows directly from the one outlined in Eqs. (|55|)-([7T|). If the 
kaon is assumed to be the intermediate meson state and the A, S, and S* are the intermediate hyperon 
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states, the strange anti-quark distribution in the nucleon is 



XSn{X, 



X 

Y i y y 



:i28) 



following Eq. ( |7TD where ry are the kaon-hyperon spin-flavor SU(6) Clebsch-Gordon coeflicients. The 
kaon distribution in the nucleon cloud, fKriv), is calculated as in Eq. (|66D with M ^ K and B Y. 
An exponential form factor, Eq. (^), was used with Ae = 1.2 GeV, corresponding to A^ = 0.94 GeV. 
The strange valence quark distribution in the kaon was assumed to be identical to the pion valence quark 
distributions. In Ref. [[I09|| the SMRS [ |113|| pion structure function was used to obtain Q^). Paiva 

et al. refer to the distribution in Eq. ( |128| ) as the intrinsic anti-strange quark distribution, after the 
intrinsic charm model of Brodsky et al. [|114|| . Their intrinsic distribution is compared to the CCFR 
data and another meson-cloud model calculation [p2| based on the GRV pion distribution functions 



[l72| in Fig. |T9|. The shape of the intrinsic distribution is valence-like, as expected in such models [|114 |. 
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Fig. 19. The intrinsic x'sn distribution calculated with the meson cloud model (solid curve). The result is 
compared to the CCFR |22] data and the model calculation of [^2| (dashed curve), labeled Ref. 3 on the figure. 
Reproduced from Ref. [109] with permission from World Scientific. 

The GRV pion-based result, labeled Ref. 3 on Fig. |T^, has the opposite behavior at low x. While both 
calculations agree with the CCFR data for x > 0.3, neither can describe the low x behavior of the data. 
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Christiansen and Magnin |[11CI|| developed their approach in the context of the valon model |j70|, described 
in Section 3.1.2. In this case, the gluon emitted by the valon can decay into an ss pair. This three valon 
+ ss state can then rearrange itself to form a kaon-hyperon bound state. They point out that the 
strange meson and baryon distributions inside the nucleon are not independent but are constrained by 
the requirements of zero strangeness in the nucleon and momentum conservation of the meson-baryon 
state, 

1 

J dx [Psix) - Pm(x)] = (129) 



1 

J dx [xPb{x) + xPMix)] = 1 . (130) 



Both conditions are satisfied if Pm{x) = Pb{^—x) where Pm{x) is the strange meson distribution in the 
proton, corresponding to Jky in Eq- ( |128|) . Note that in this model, no explicit form factor is required 
since the valon model is defined on the light cone. Thus, the conditions in Eqs. ( p.29| ) and ( |13CI| ) can be 
satisfied. 

They use the recombination model |]T3[ in Eq. (^) with a = 1 in Eq. ( |75D to predict the rearrangement 
of the valons and the strange partons in the proton so that 

Pm{x) = f -F{y,z)R{x,y,z) . (131) 

J y J z 

The normalizations of F{y, z) and R{x, y, z) are fixed by assuming the probability of ss pair production 
by a valon is ~ 4 - 10% pO| . 

The nonperturbative strange and anti-strange quark densities in the nucleon are 

^^^(x) = / ^Pniy) SB (-] (132) 
5^^(x)=/^PM(y)^A/f-') . (133) 

i y \yJ 

Two different approximations of sb{x) and sm{x) were tested, 

sm{x)=Qx{1-x) sb{x) = l2x{l -xf (134) 

sm(x) =g"(a;)(l - xf-'' sb{x) = ^u;ix) . (135) 
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The strange and anti-strange densities in the proton, Eqs. ( |132| ) and (|133|) , are shown in Fig. ^ along 
with the difference — Ahhough they do not compare the calculations directly to data, it 

is clear that the kaon and hyperon distributions in Eq. ( p.35| ) agree better with the trends of the CCFR 
data since the differences between the s and s distributions are small except at low x where the valon 
model is expected to break down. An "intrinsic" -type meson cloud model calculation was performed 
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Fig. 20. The upper plot shows the nonperturbative strange (thin lines) and anti-strange (thick lines) distri- 
butions in the nucleon. The solid curves are calculated assuming the schematic kaon and hyperon distribu- 
tions of Eq. ( |134| ) while the dotted curves are calculated with Eq. ( |135| ). The lower plot gives the difference 
s^^ {x) — s^^ {x) obtained with Eq. (|134| ) (solid curve) and Eq. ( |135|) (dotted curve). Reproduced from Ref. [110] 
with permission from Elsevier Science. 



by Brodsky and Ma ||111|| with results similar to those shown in the solid curve in Fig. |20|. A light-cone 
calculation of s{x) — s{x) by Melnitchouk and Malheiro shows the opposite behavior: s{x) — s{x) is 



positive at low x and negative at large x [|112|| 
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Ji and Tang discussed the implications of the CCFR data on the "locahty" of the strange sea ||107| 



They define two hmits on the locahty. If the s and s are bound in pairs or if they move independently 
of each other but have similar interactions with other partons, the s and s have essentially the same 
spatial and momentum wavefunctions. This limit corresponds to small locality. If, instead, the quarks 
and anti-quarks move independently of each other with different interactions with other partons, their 
distributions should be different, corresponding to large locality. The locality in coordinate space, the 
strangeness radius, is 

(r^) = / dfr^ {\s{f)\^ - |s(f)n . (136) 



The corresponding locality in momentum space is |p,07 
1 

dx \s{x) — s{x)\ . (137) 







A small Lg or (r^) is equivalent to a large locality while a large Lg or (r^) represents a small locality. 
They concluded that the standard meson cloud picture of the s and s distributions, as in Eq. (I 



|62| , |109|] , predicts a larger locality than the CCFR data can accommodate, see Fig. ^ They attempted 



to model the s and s distributions in a manner that could correlate coordinate and momentum space 
localities. One obvious way of doing this is to assume that g —>■ ss is the dominant ss production 
mechanism at scale /x^ = 1 GeV^. This perturbative production is regulated by introducing effective s 
and s masses which could arise from mean-field interactions between the s and s with other components 
of the sea. They allow ^ rris to account for differences between the s and s interactions in the nucleon. 
The larger the effective mass, the harder the momentum distributions. They calculate the strange quark 
distribution in the gluon by |p.07| 



fs{x,^i^) = ^J dkl[{kl + m|) + (1 - x){ml - m|)]-' {2x{l - x)msms 



+ {kl + m|)[x2 + {x- If] + (1 - xf{ml - m|)} (138) 
where as = 0.5 and f-s{x,fi'^) = /^(l — x, /i^). 

The strange quark distribution in the nucleon is obtained by convoluting /^(x, /i^) with the nucleon 
gluon distribution, 

— xfs{-, nis, mi) G{y, = 1 GeV^) . (139) 

. y y 



54 



They fit Eq. (|139|) and the corresponding anti-strange distribution to the CCFR s and s results by 
adjusting and m^. With the CTEQ3 gluon distribution | ]115| , they found = 260 ± 70 MeV and 
TiTs = 220 ± 70 MeV. Note that the lower limit of rris is correlated with the upper limit of rris- Within 
the range of uncertainties, the effective masses are equivalent and the agreement with the data is good, 



as seen in Fig. 21 
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Fig. 21. Comparison of the strange and anti-strange quark distributions calculated with Eq. (|139|) with the 



CCFR data [^2[. The effective masses = 260 MeV and = 220 MeV are used. Reproduced from Ref. [107] 
with permission from Elsevier Science. 



Comparing the meson-cloud type models of Figs. |T9| and [2^ with the gluon splitting picture of Fig. 
it appears that the s and s distributions are more strongly correlated than allowed by the meson cloud 
model unless rather hard valence distributions are assumed for the valence quark distributions in the 



strange mesons and baryons, such as in the dotted curves of Fig. ^ The data thus seem to imply that 
the s and s quarks undergo similar interactions and their distributions have a high degree of locality. 

The rather significant correlation between the s and s distributions suggests that the strangeness radius, 
Eq. ( |136| ), is small. The strangeness radius can be extracted by measuring the strange quark contribution 
to the nucleon elastic form factor, as done by SAMPLE |9^. Both Paiva et al. [ |109|| and Ji and Tang 



||107|| provided estimates of the strangeness radius within their models. 



In the intrinsic strangeness calculation of Paiva et al. the integral over sat is equivalent to the probability 
of finding an intrinsic ss pair in the proton wavefunction. They obtain Pis = / dxs]\f{x) = 0.12. The 



55 



probability is related to the strangeness radius by 
\{r!)\ 



IS 



(140) 



where, using = 0.72 fm and the integrated value of Pis, they find ||109 
\{rl) \ = 0.0622 W . 



(141) 



They also note that if they only considered KA components in the cloud, the strangeness radius would 
be reduced by more than a factor of two. However, they can only predict the absolute value of | (r 
not its sign. 



Ji and Tang translate their small locality in momentum space into a coordinate space locality by 
dimensional analysis ||107 



1(^^)1 



< 0.005 fm^ 



(142) 



considerably smaller than the prediction of Paiva et al. [|109|| as well as other earlier theoretical predic- 
tions [p.04]J107|| although the result is similar to that of Ref. [|112|| . 



3. 2. 3 Strange form factors 



The strange form factors have been studied phenomenologically in the context of the E734 u and u 



scattering experiment |Q. The neutral current cross section is sensitive to the electric, magnetic, and 
axial strange form factors of the nucleon while the charged current interactions do not access the strange 
quarks in the nucleon sea at low Q^. 



The neutral current elastic scattering cross sections are [pJi 
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(144) 
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where ew = ^ — 4sin^6'vy and Gp is the Fermi constant. The dependence of the form factors has 
been suppressed in Eqs. ( [I43| ) and ( |144D . The kinematic quantities Cm, Ce, Ca, and Cam above are 
defined as 



CM = ly' 
Ce = 1 
Ca = 1 



y 



My 
1 



Cam = 2y(l- -y 



^ 2^ My 
1 



(145) 
(146) 
(147) 
(148) 



where y = Q^/ {2mpE^) and E^, is the incident neutrino energy. The neutral current cross sections include 
terms independent of the strange form factors as well as terms with linear and quadratic dependencies 



on Gli, G%, and F|. 



The corresponding charged current cross sections are [116 



Or p 



„ , Vud\ \ '^M Xy'M ~ ^M) + '-'E ['^E ~ '^E ) 

+ CaFI + CamFa{GZ-G'I:})} 



— — — — \^ud\ \'^M\S'm~^m) +tv-El<-f£ — tj 



dQ" 



771 \ 2 
E I 



+ CaF\ — CamFa {GYi — C]/ )| 



(149) 



(150) 



where |VudP is the square of the Cabibbo-Kobayashi-Maskawa matrix element mixing up and down 



quarks. The dependence of the form factors has again been suppressed in Eqs. ( |149| ) and ( |150| ). It is 
necessary to check whether the equivalent free scattering cross section obtained by E734 for quasielastic 
knockout is appropriate to use to extract the relative strengths of the strange magnetic form factors 
from the E734 charged current data. Alberico et al. showed that while including final-state interactions 
between the ejected nucleon and the residual nucleus suggests that only half of the ^^C reactions 
correspond to elastic proton knockout, the ratios, Eqs. ( |125| )- ([T27|) are affected at the percent level for 
neutrino energies greater than 1 GeV ||117|| . Since E734 reported ratios, Eqs. ( |125|) - (|127|) , only elastic 



scattering of neutrinos on free protons were considered. 

Integrating the cross sections over Q^, it is possible to compare these results with the ratios reported 
by E734. The neutral current ratio R, Eq. ( |125|) , is formed by the ratio of Eq. ( |144|) to Eq. ( |143| ) . The 
charged current ratios R^, Eq. ( |126D , and Rtj, Eq. ( |127[ ), are the ratios of Eq. ( |143[ ) to Eq. ( |149| ) and 
Eq. dm) to Eq. respectively. 
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By forming the asymmetry of neutral current to charged current interactions, defined as the ratio of 
the difference between Eqs. ( |143[ ) and ( |144D to the difference between Eqs. ( p.49| ) and 



{dayp^yp/dQ'^) - {da-pp^yp/dQ'^ 



(151) 



(dcTyn^^-p/dQ^) - {dayp^f,+n/dQ^) 
model independent information on the axial, F^, and magnetic, G^, strange form factors of the nucleon 



can be obtained 117 . The differences in the cross sections can be written in terms of the form factors 



so that the asymmetry becomes 
1 



2 sin^ 9w 



M 



(152) 



where Fa is the charged current axial form factor and 



0.5(Gf 



— ^J.u<^'^JJ^J — G]^) is the charged current 
isovector magnetic form factor of the nucleon. Note that A is independent of the electric form factors. 

To compare the calculations to the neutrino energy and momentum transfer averaged BNL E734 ratios, 
in Eqs. ( |125| )- ([T27|) , they integrated over the neutrino energy spectrum, Eq. ( |124| ), and to obtain the 
integral asymmetry 



(A) 



(.(^ vn~* ii~ p) {^T/p^/i+n} 

Rujl-R) 

1 RRy I RjJ 

0.136 ±0.008 ±0.019 



(153) 

(154) 
(155) 



from the data, Eqs. (|T25|)-(p7|). 

Dipole parameterizations of the non-strange and strange form factors were assumed. The strange axial 



and vector form factors are 11117 



(Ml ± Q2)2 

m + Q'f 
Ami (M^ + Q^y 



(156) 
(157) 
(158) 



where the strange couplings g%, fig and pg are left as free parameters. The vector and axial cutoff masses 
are My = 0.84 GeV and Ma = 1.032 GeV respectively. Other parameterizations of the form factors 
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stronger than the dipole would reduce the overall effect of strangeness at high but in the range 
0.5 < < 1 GeV^ changing the dependence has a very weak effect on the results. 



The ratios Ry^ Rtj, and R and the asymmetry A are shown in Fig. ^ along with several choices of g%, 
/is, and ps- Note that A is independent of pg. The ratios and asymmetry are given as a function of ps for 
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Fig. 22. The ratios R^, Rjj (a), and R and the asymmetry A (b) are calculated as a function of /Ug and compared 
to the E734 data (shaded areas). All the calculations correspond to elastic scattering. Two choices of the axial 
coupling g\ are shown: and —0.15. For each value of g\, three values of ps are chosen: —2 (dot-dashed), 
(solid), and 2 (dashed). Reproduced from Ref. |117] with permission from Elsevier Science. 



g\ = Q and —0.15 with ps = —2, 0, and 2. Changing g\ influences the ratios more than does changing 
Ps, especially for Ry. In fact, the sensitivity of R-jj to g% excludes g\ < —0.15. The calculations of R^ 
and A are most consistent with the E734 data when g% is close to zero and ps is negative. Varying the 
axial cutoff Ma by ±3.5% only has a significant effect on R-jj. Lowering Ma increases the allowed range 
of g\ if /is = Ps = while increasing Ma restricts g\ to values close to zero. 

Although the results obtained by the SAMPLE collaboration, Eq. ( |121| ), suggest that ps should be 
positive, the large uncertainties cannot rule out Ps < 0. Thus so far the E734, SAMPLE, and HAPPEX 
results are consistent with each other. As pointed out in Ref. [ 117|| , if the parity-odd asymmetry mea- 
surements are improved, then the more stringent bounds on ps can place better constraints on g\ and 
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ps- The strange form factors have also been computed in a hght-cone approach ||118|| and compared to 
the HAPPEX data. In this case, Melnitchouk and Malheiro calculated the combination of G% and G%i 
in Eq. (|123| ) and found that the strange electric and magnetic form factors were small and positive, 
consistent with the experiments. 



3.2.4 Strange sea summary 



The difference between the s and s distributions determined by the CCFR collaboration has been 
studied in the context of meson cloud models [ p.09| , |110| , [lll| , |112|| and has been modeled by gluon splitting 



into s and s quarks with different effective masses |107|. Meson cloud models tend to produce larger 



differences between s{x) and s{x) than can be accommodated by the data, leading Ji and Tang to suggest 



that the strange sea is highly localized [107]. However, further measurements are needed with analyses 
based on more recent parton distribution functions. A reanalysis of the neutrino data is underway which 



could provide some answers |95 



Information on strange form factors in the nucleon has been obtained by the SAMPLE [PTl , HAPPEX 



98 1, and E734 |99| collaborations. The results are consistent with each other and show that the strange 



form factors could be small and positive but are also consistent with zero [117 . 



3.3 Heavy quark contributions to the sea 



Quarks more massive than the strange quark, e.g. charm and bottom quarks, have also been included 
in global analysis of parton distribution functions. (See Ref. [^ for a recent discussion and further 
references.) Early treatments of heavy quark parton distributions assumed that below a scale ~ mg, 
there are rif flavors in the perturbative sea and in the calculation of a^. When > ttIq, rif + l massless 
flavors are included in the perturbative sea and the running coupling constant. Thus above 
the heavy quark is treated as massless 



m 



Q' 



2q] . This scheme is now generally referred to as the Variable 



Flavor Number scheme. At the opposite extreme, the heavy quark is never treated as part of the nucleon 
sea but is produced perturbatively through photon-gluon fusion. The number of flavors remains flxed, 
regardless of Q^. This treatment, used in the GRV 94 parton densities [27|, is referred to as the Fixed 
Flavor Number scheme. Obviously neither of these schemes is correct at all Q^. The Fixed Flavor 
Number scheme is good near threshold but cannot incorporate large logarithms at ^ mq while the 
heavy quark should not be treated as massless near ~ mg, as in the Variable Flavor Number scheme. 
Interpolating schemes, such as the one used by CTEQ []25| , |119|| , have been introduced which reproduce 
relevant features of both schemes. The heavy quarks are essentially as produced by photon-gluon fusion 
when ~ mg and behave as massless quarks when S> rriq. All three schemes have been included 
in the recent global analysis by the CTEQ collaboration. The charm quark distribution from three of 
their flts is shown in Fig. E3|. 
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Fig. 23. Charm quark parton distribution in the Variable Flavor Number scheme, CTEQ 5M (solid), the 
Fixed Flavor Number scheme, CTEQ 5FF (dashed), and the interpolating scheme of Ref. [[119(| , CTEQ 5HQ 
(dot-dashed). 



Charm DIS data are discussed first. The EMC large x charm data ||120|| imphed some nonperturbative 
source of charm production. These data have not been followed up with subsequent large x charm 



measurements. More recent HERA data ||121|,|122|| suggest that only perturbative production is needed 
at low X. Models of nonperturbative charm production are then introduced which could be important 
at large x. 



3.3.1 Current experimental evidence 

The European Muon Collaboration, EMC, [|120|| at CERN first measured the charm structure function. 



F2c- In the relatively large x region of the EMC measurement, charm only contributed ~ 1% to the 
total structure function F2^. Charm was identified by dimuon events, where one muon was consistent 
with charm meson decay, and trimuon events, which were consistent with a pair produced by 

the semileptonic decays of a pair of charm mesons. In both cases, the additional final-state muon was 
the scattered projectile muon. They compared their results to a leading order photon-gluon fusion 
calculation. In their analysis, they assumed a simple gluon momentum distribution, xg{x) = 3{1 — x)^ 
and used the cc pair mass in the scale of the strong coupling constant, fi"^ = Q'^ + m^^. They included 
~ 83% of the cc cross section below the DD threshold in the total charm cross section. Their data, 
which has never been incorporated into global analyses of parton distributions, are shown in Fig. |2^ as 
a function of with uncertainties on both x and given. A nonnegligible dependence is clearly 
indicated. 

More recently, the charm structure function has been measured to very low x at HERA ||121| , |122|| . At 
these low x values, Fg^f/F^^ ~ 0.25 at ~ 11 GeV^, a surprisingly large fraction. The charm structure 
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Fig. 24. The charm structure function as a function of measured by the European Muon Cohaboration 



pO| |. From top to bottom the values of the x bins are: 0.0044, 0.0078, 0.014, 0.025, 0.044, 0.078, 0.14, 0.25, and 
0.44. To separate the data, Fg^^ in each x bin is scaled by a factor of 4 from the next higher x bin. Therefore, 
only the highest x data (lowest points) have the correct scale. Only the statistical uncertainty is shown. 

function is measured in the reactions e~^p — > e^D*^X by studying the decay channels D*~^{2010) 
Z}°(1864)7r^ + charge conjugate where 7r+ is a slow 7r+. The and from D*^ decays subsequently 
decay into K{n7c) final states, —>■ K^tt^ or K^ti^ti^ti^ (+ charge conjugates) resulting in D*^ 
decays with two or four pions in the final state. Once the final states have been reconstructed, the cc 
cross section is obtained from the D*^ production cross section, extrapolating to the full kinematical 
phase space, and employing the hadronization fraction F{c D*^) to obtain the total charm cross 
section. The extrapolations necessary to go to full phase space were significant, resulting in a factor of 
ten correction for the four-pion final state at low and a factor of four at high Q^, and neglect charm 
production by processes other than photon-gluon fusion. The cc production cross section is then 



g2) 27va^ 



{ [l + (1 - yy] Q') - y^Fllix, Q^)} 



dxdQ'^ xQ'^ 
where the contribution from the second term is negligible. Equation 



(159) 



can be obtained from Eq. (|TTp 



for the virtual photon by changing integration variables from y to Q and taking the rrih/E <^ 1 limit. 
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The hadronization fraction of charm to D*, assuming that final state production of charm hadrons is the 
same in e+e" annihilation and DIS, is F{c D*) = 0.222 ± 0.014 ± 0.014 |T5B[. They assume = 1.4 
GeV and the renormalization and factorization scales used in the photon-gluon fusion production of 
charm at NLO are = ^4771^ + Q'^. Recall that the EMC data were analyzed with fi = ^ m^^ + Q"^ . 
The data from the two different decay modes are combined in the final data set and the charm structure 
function is extracted from Eq. ( |159| ). The ZEUS data are shown as a function of in Fig. |25|. The 
structure function shows significant dependence as well as a steep rise as x decreases with 
fixed. The ZEUS data and the low x EMC data are consistent with each other and with perturbative 
generation of the charm sea. 



Since HI has a smaller charm sample, they only present the ratio F2^/F2^ averaged over the kinematic 
range 10 < < iqo GeV^ and 0.0008 <x < 0.008. They found F:JI/ F:J^ = 0.237 ± 0.021 toml an 
order of magnitude greater than at larger x. Their result is compatible with the ZEUS result in the 
same kinematic interval. 
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Fig. 25. The charm structure function as a function of measured by ZEUS |121|. From top to bottom 
the values of the x bins are: 0.00005, 0.00013, 0.0003, 0.0005, 0.0008, 0.0012, 0.002, 0.004, 0.008, and 0.02. To 
separate the data, in each x bin is scaled by a factor of 4 from the next higher x bin. Therefore, only the 
highest X data (lowest points) have the correct scale. Only the statistical uncertainty is shown. 
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3.3.2 Nonperturbative heavy quarks 



The inclusion of heavy quarks in the perturbative generation of the sea and global analyses of the parton 
distributions has already been discussed. However, the large x and EMC charm structure function 
data hinted that these results were inconsistent with production of charm by photon-gluon fusion. This 
led to the speculation that charm could be generated in the sea by a higher-twist mechanism. Two 
nonperturbative models of charm in the sea are discussed here, intrinsic charm | p. 14] ] and the meson 
cloud model |09l|l2l. 



3.3.3 Intrinsic charm 

The QCD wavefunction of a hadron can be represented as a superposition of quark and gluon Fock 
states. For example, at fixed light-cone time, a hadron wavefunction can be expanded as a sum over 
the complete basis of free quark and gluon states: |\E'/j) = J2m \^) i^m/h{xi, kx^i) where the color-singlet 
states, |m), represent the fluctuations in the hadron wavefunction with the Fock components \qiq2(l3), 
\QiQ2<l3g)i |9i?293cc), etc. The boost-invariant light-cone wavef unctions, 4'm/h{xi, kT,i) are functions of the 
relative momentum coordinates Xi = kf / and kT,i- Momentum conservation demands Z^iLi ^« — 1 
and J2i'=i^T,i = where n is the number of partons in state |m). For example, intrinsic charm, IC, 
fluctuations ||114| | can be liberated by a soft interaction which breaks the coherence of the Fock state 



provided the system is probed during the characteristic time that such fluctuations exist. 



Microscopically, the intrinsic heavy quark Fock component in the proton wavefunction, \uudcc), is 
generated by virtual interactions such as gg —>■ QQ where the gluons couple to two or more valence 
quarks. The probability for cc fluctuations to exist in a hadron is higher twist since it scales as 
relative to the extrinsic, EC, leading-twist production by photon-gluon fusion | |12(i| . 



The dominant Fock state configurations are not far off shell and thus have minimal invariant mass, 
M'^ = J^iTri^/xi where = mf + (/cf^j) is the square of the average transverse mass of parton i. The 
general form of the Fock state wavefunction appropriate to any frame at fixed light-cone time is 

^(^^-M = ^^^ (160) 



mt 



where F is a vertex function, expected to be a slowly- varying, decreasing function of m\ — M^. The 
particle distributions are then controlled by the light-cone energy denominator and phase space. This 
form for the higher Fock components is applicable to an arbitrary number of light and heavy partons. 
Intrinsic cc Fock components with minimum invariant mass correspond to configurations with equal 
rapidity constituents. Thus, unlike extrinsic heavy quarks generated from a single parton, intrinsic 



heavy quarks carry a larger fraction of the parent momentum than the light quarks in the state ||114 |. 
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The parton distributions reflect the underlying shape of the Fock state wavefunction. Assuming it is 
sufficient to use (k^) for the transverse momentum, the probability distribution as a function of a; in a 
general n-particle intrinsic cc Fock state is 



ic 



(161) 



where A^„ normalizes the n-particle Fock state probability. 



The intrinsic charm structure function at leading order, LO, and next-to-leading order, NLO, is now 
described. At LO in the heavy quark limit, nic, rric ^ rrih, friq, 
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(162) 



If n = 5, the minimal intrinsic charm Fock state of the proton, \uudcc)^ is assumed and Pic = 1% ||114 
then = 36. In what follows, a 1% intrinsic charm contribution is assumed unless otherwise noted. 



Hoffmann and Moore ||127|| incorporated mass effects into the analysis. They first included threshold 
effects by introducing a scaling variable ^ = 2ax/[l + (l+4x^mp/(5^)^/^] with a = [l + (l+4m^/Q^)^/^]/2. 
The cc mass threshold requires ^ < 7 < 1 where 7 = 2ax[l + (1 + Ax'^rrip/ Q'^Y^'^]^^ and x = Q'^/[Q'^ + 
Ami ~ ^p] ■ Then Eq. ( |162D is replaced by 



^,^^^°(x,Q^m^) = ^ec(e,7) 



(163) 



with c(2;,7) = c{z) — zc{'y)/'y for 2; < 7 and zero otherwise. The final LO result can be written more 
generally as in Eq. (18) of Ref. [ |127|| , 



F,^?^0(x,Q^m; 



2 _2\ 



8x2g3 



9(g2 + 4x2m2)3/2 



^(^' ^) + f (^' ^) 



(164) 



where 

^(e,7) 



Gxm^ 



{Q^ + 4x^ml)J t 



dt , , / m 
■c(t,7) 1 



2xtml 2xml 
1 + + 



The NLO IC component of the structure function is given by 
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pic NLO , 



€/7 



(165) 



The lowest order cross section is normalized to a^2\^i'^) 



6{1 — z). The NLO QCD corrections to 
the IC contribution are given in Eq. (51) of Ref. ||127|| . The IC results at NLO are the sum of the 
kinematically corrected LO result, Eq. ( [164| ), and the full NLO correction. 
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Fig. 26. (a) The contributions to the structure function F2^{x,Q'^,m'^) at = 7 GeV^: the massless result, 
Eq. ( [L62D , (upper dotted line); the ^-scaling result, Eq. ( |163| ), (dot-dashed line); and the kinematically cor- 
rected formula, Eq. (|164D , (top solid line). Also shown are the NLO corrections given by Eq. ( |165D with the 



leading-log result, Eq. (54) i n [|127[| , (dashed line), and the full result, Eq. (51) in |127], (lower dotted line). 
The sum of Eqs. (|l6|) and {^6^ using the full result, Eq. (51) in ||127| ], represents the total IC contribution 
to F^^^^'^{x,Q'^,ml) (lower solid line), (b) Same as part (a) for = jq GeV^. Reproduced from Ref. |2| 
with permission from Elsevier Science. 



7 and 70 GeV^. At LO the 
), and the full calculation, 



Figure ^ shows the LO and NLO IC contributions to F2{x, Q , ml) for Q 
massless, independent, result, Eq. ( |162| ), the ^ scaling formula, Eq. (|16 
Eq. ( |164| ), show substantial differences at low while the three are nearly indistinguishable when 
is large. The NLO corrections calculated using Eq. ( |165[ ) are also shown. The leading-logarithmic 
approximation, Eq. (54) in ||127|| , and the complete result, Eq. (51) in |[127|] , differ at small but 
are almost identical at large Q^. The total NLO IC results are actually negative at large x and Q^, 



indicating that even higher order terms are needed in the perturbative expansion. 
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The EMC data have been analyzed with a combination of photon-gluon fusion and intrinsic charm, first 
by Hoffmann and Moore [|127|| and later by Harris et al. [ |128| |. In both cases, the intrinsic charm calcu- 
lation was carried to next-to-leading order. While only leading order photon-gluon fusion was available 

in the more recent analysis, the photon-gluon fusion calculation was carried 
Photon-gluon fusion is used to generate charm perturbatively 



for the first study ||127 



128 



to next-to-leading order as well 
because at the measured of the EMC experiment, the mass of the charm quark is not small. 



To NLO ||129|| , the extrinsic charm structure function is calculated from the inclusive virtual photon- 
induced reaction at the parton level, 7* -|- ai ^ c + c + a2 where Oi and 02 are massless partons. In an 
inclusive calculation, the c and 02 momenta are integrated over. The charm structure function is the 
sum of the longitudinal and transverse virtual photon total cross sections, oc {Q"^ /4:7r^a){aT + o"l). 
At NLO, the extrinsic charm structure function is then 



«0 



Co 



i=q,q 



%9 



"'2,j "T U2^i 



(166) 



where ^0 = a;(4m^ -|- The quark charges are in units of e. The parton momentum distributions 

in the proton are denoted by where /i, the mass factorization scale, has been set equal to the 

renormalization scale in a^. The scale independent parton coefficient functions, i, / 



O5 ^2 ] , (il] 



and 0^2} originate from the coupling of the virtual photon to the partons. The functions and 0^2] 
represent the virtual photon-charm quark coupling, thus appearing for both charged parton and gluon- 
induced reactions. The virtual photon-light quark coupling gives rise to d^2 \ while 02] is an interference 
term proportional to CcCj. 



The EMC data ||130|| with V = Q"^ /2mpX = 53, 95, and 168 GeV were fit by a sum of the photon- 
gluon fusion and IC components ||128|| . The normalization of both the IC and EC components are free 
parameters so that 



(167) 



with the scale /i = yml-^ + Q"^. Since /i(m^^), an exclusive cc production calculation at NLO | [L31|| was 
carried out to more precisely match the EMC analysis. The parameter e, typically larger than unity, can 
be considered an estimate of the NNLO contribution to photon-gluon fusion. Since a 1% normalization of 
the IC component is already assumed, the fitted 5 is the fraction of this normalization. Several different 
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u = 53 GeV 



u = 95 GeV 



u = 168 GeV 



PDF 



CTEQ3 1TT|] 1.0 ± 0.6 0.4 ± 0.6 1.2 ± 0.1 0.4 ± 0.3 1.3 ± 0.1 0.9 ± 0.5 
MRS G [|l3gi 1.0 ± 0.7 0.3 ± 0.6 1.4 ± 0.2 0.3 ± 0.3 1.5 ± 0.1 0.8 ± 0.5 
GRV 94 pi 1.2 ± 0.8 0.3 ± 0.6 1.4 ± 0.2 0.3 ± 0.3 1.5 ± 0.1 0.9 ± 0.5 



Table 4 

Results of the least squares fit of the photon-gluon fusion and IC contributions to the EMC data according to 



Eq. ( 11671) . Adapted from p28 |. 



parton distribution functions were used which produced essentially the same results. The results are 
presented in Table ^. The errors in the table correspond to a fit at the 95% confidence level. The final 
results for the combined model of Eq. ( p.67| ) are shown in Fig. for the CTEQ3, MRS G, and GRV 
94 parton densities. Given the quality of the data, no statement can be made about the intrinsic charm 
content when z/ = 53 and 95 GeV. However, with z/ = 168 GeV an intrinsic charm contribution of 
(0.86 ± 0.60)% is indicated. These results are consistent with those of the original analysis by Hoffman 
and Moore P7| . 
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Fig. 27. (a) The EMC data for the structure function F2c{x, Q'^,mf.) at u = 53 GeV are plotted as a function of 
X together with the fit from Eq. (|l6|) with e and 6 taken from Table g The CTEQ3 (sohd), MRS G (dashed) 
and GRV 94 HO (dotted) parton densities are shown, (b) Same as (a) for = 95 GeV. (c) Same as (a) for 



ly = 168 GeV. Reproduced from Ref. [128] with permission from Elsevier Science. 



68 



Only the intrinsic charm contribution to the charm structure function has been discussed here. There are 
many other apphcations of intrinsic charm in charm hadron production. See e.g. | |1 26| , |1 3^j| , p, 34| jl 35| , |1 36| | 
for more details. 



3.3.4 Meson cloud models 

Paiva et al. have also calculated an "intrinsic" charm component of the nucleon sea within the context of 



the meson cloud model | 109 |. They assumed that the nucleon can fluctuate into DA^. The c distribution 
in the nucleon is then 



1 



XCn{X 



)= [dyfM^-OnO- (168) 



X 



where 



2 ^max r / 



\21 



—00 ^ 

and tmax = mj^y — m\^y/{l — y). In this case they chose a monopole form factor with = 1.2 GeV. 



The coupling constant was assumed to be gjjj^f^^ = —3.795. From heavy quark effective theories ||137 
the c distribution in the D is expected to be hard because in the bound state, the c exchanges momenta 
much less than m^. They make the extreme assumption that the entire D momenta is carried by the 
charm quark, cjy = x6{x — y). The resulting c distribution is compared to the intrinsic charm quark 
distribution in the heavy quark limit, Eq. (|162|) , with xcn = xcn = 9F2^^'^/8 in Fig. ^ 



More recently Steffens et al. investigated all the charm structure function data with two variants of 
intrinsic charm ||124|| . The first was that of Eq. ( |162| ), called ICl in their paper, while the second was 



a meson cloud model, IC2. In the second approach, the c distribution is obtained from the light cone 
distribution of D mesons in the nucleon. 



c'^'^x) ^ fjj{x) (170) 
1 7^^2 a^i^^ kl) k\ + (mA, - (1 - x)mM f 



IGTT^y x{l- x){s-^^^-m%f l-x 
A hard charm momentum distribution was assumed in the similar to that of Ref. |109|. The vertex 



function k\) is parameterized as g"^ = g^^A^ + m^)/(A^ + s-^j^J where 5-5^^ is the square of the 
center of mass energy of the -DAc system and g^ the coupling constant at sj^^^ = m^. For an intrinsic 
charm probability of 1%, A ^ 2.2 GeV. The charm distribution is then 
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Fig. 28. The cat distribution calculated with the meson cloud model, Eq. (|168D , (solid curve) compared to 
the intrinsic cat distribution of Eq. (|162| ) (dashed curve) labeled Ref. 19. Reproduced from Ref. [109| with 
permission from World Scientific. 



3 /3a;\ 



'1711 



where the charm distribution in the Ac is assumed to be ~ S{x — 2/3) and /a^(x) = fjy^l — x). 



Using several different schemes to calculate the extrinsic charm sea, Steffens et al. compared their two 
component model to the combined ZEUS and EMC data. As expected from the agreement of the low 
X ZEUS data to the perturbative generation of the charm sea, the ZEUS data show no sign of intrinsic 
charm. However, the EMC data still indicate a need for intrinsic charm at large x and Q^, as found in 



Ref. ||128|| . The best results, shown in Fig. |2^, are obtained in the fixed flavor scheme, as outlined at the 
beginning of this section. They find Pic ~ 0.75% for the standard massless intrinsic charm at leading 
order, their ICl, and Pic ~ 0.4% with their meson cloud approach, IC2 ||124|| . They claim that the data 
indicate a slight preference for model IC2. However, unlike the results shown in the previous section, 
they did not include the NLO corrections to which reduce the large x structure function relative 
to Eq. (|162|) , as seen in Fig. |2^. Clearly, more data are needed at large x and to finally resolve the 
situation. 
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Fig. 29. The charm structure function calculated in the fixed flavor number scheme (FFNS) for three bins 
with the intrinsic charm models ICl, Eq. ( [162D , and IC2, Eqs. ( |17[1| ) and ( |171| ). The ICl result is normalized to 
0.75% while IC2 is normalized to 0.4%. The most recent GRV parton distributions |138|] are used for the light 



quark and gluon densities. Reproduced from Ref. |124] with permission from Springer- Verlag. 
3.3.5 Heavy sea summary 

Considerable progress has been made in the understanding of the charm component in the nucleon sea, 
particularly at small x. Treating the heavy quark as massless at large and massive at ~ uiq put 



the heavy quark contribution to the sea on a firmer theoretical footing. Charm data at HERA |121,122 
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are in good agreement with perturbative production of charm. However, the large x and EMC data 
| ]120 still do not fit into this scheme. Intrinsic charm on the level of ~ 1% can explain this lower energy 



data 1 114, 128 1 . Whether the c and c distributions are identical |114| or differ as in an e.g. DAc cloud 



||124|| is unclear. More data in this interesting region are needed. 



In addition to the reconstructed charm structure function data at HERA, deep- inelastic scattering 
data with significant statistics at high and high x are now available ||139| , |14CI|| which could also 
shed some light on charm production processes. There is an excess in the data beyond the predictions 
based on standard QCD parton distribution functions. These data have been compared to calculations 
with intrinsic charm that enhance the standard model charm cross section ||141| , |142|| in the x and 
regions of interest to the HERA excess. It was found that intrinsic charm could account for a significant 
enhancement over the standard model predictions in charged current scattering. More data are needed 
to confirm the effect and to clarify the role of charm in this region. 

A nonperturbative bottom sea component in the nucleon has not been discussed here but would be a 



plausible extension of the models. See Ref. [126| for some predictions of intrinsic bottom effects. 



4 Summary 



There is only room in this review for discussion of a limited subset of all experimental results and 
theoretical progress towards understanding the data. For more in-depth discussion of the work described 
here as well as other results, the reader is encouraged to consult the appropriate references. 

As discussed in this review, studies of the nucleon sea have revealed a rich structure compared to early 
evaluations of parton distribution functions with an SU(3) symmetric sea, e.g. Ref. I^Sj. Of these studies. 



the Gottfried sum rule violation, which showed that ci > m, is the best documented experimentally. 
Interpretation of the NMC and E866 data as an indication of a cloud of virtual mesons and baryons 
surrounding the nucleon has proved fruitful. 

Fewer data are available so far for the more massive components of the sea such as strangeness and 
charm. While the CCFR data suggest that a difference between the s and s distributions cannot 
be excluded, the distributions appear to be more strongly correlated than a kaon-hyperon cloud can 
accommodate. The high x and EMC data still indicate the presence of a nonperturbative intrinsic 
charm component of the nucleon sea but the exact form of such a component is not fixed. In both cases, 
more data are needed to resolve both the origin and shape of the strange and charm contributions to 
the sea, particularly at large x and Q^. 

A relatively short time ago, there was a strong bias that the nucleon sea was flavor symmetric and that 
the quark and anti-quark distributions should be identical. These expectations have been exploded as 
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experiments are able to probe individual quark and anti-quark distributions more precisely. Increasingly 
good data and more sophisticated analysis techniques can hopefully resolve some of the issues presented 
here. This field is obviously very data-driven. More surprises undoubtedly await those who dive deeper 
into the "cloudy" depths of the nucleon sea. 
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